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ABSTRACT 


Capillary Flow in an Interior Corner 


Mark M. Weislogel 


The design of fluids management processes in the low-gravity environment of space 
requires an accurate model and description of capillarity-controlled flow in containers 
of irregular geometry. Here we consider the redistribution of fluid along an interior 
corner of a container; following the rapid reduction of gravity, —gz, the fluid advances 
along the corner to a height z — The analytical portion of the work presents an 
asymptotic formulation in the limit of a slender fluid column, slight surface curvature 
along z, small inertia, and low gravity. The scaling introduced renders all parameters 
in closed form or containing a slightly varying coefficient which may be computed 
numerically. This scaling clarifies the distinction between the effects of corner ge- 
ometry and surface curvature. Introducing an invariant transformation allows the 
asymptotic equations to be solved using an accurate and efficient numerical scheme. 
New similarity solutions are found and a list of closed form expressions is provided 
for flow rate and column length, C. In particular, it is found that C is proportional 
to t 1 / 2 for the constant height boundary condition, t 2//5 for a spreading drop, and t 3 / 5 
for constant flow. A linear solution which describes the tip of the fluid column is 
found to exist for all solutions. In the experimental portion of the work, measure- 
ments from a 2.2s drop tower are reported. An extensive data set, collected over a 
previously unexplored range of flow parameters, includes estimates of repeatability 
and accuracy, the role of inertia and column slenderness, and the effects of corner 
angle, container geometry, and fluid properties. It is found that there exists a loca- 
tion zh at which the interface height remains constant h(zn,t) = H after an initial 
transient, t>t H - The height H is equal to the value predicted from the static shape 
of a constant curvature fluid interface. Comprehensive comparison is made between 
the analysis and measurements using the constant height boundary condition. Good 
agreement is found for all times t > tn except near the tip for which agreement is 
achieved only as t — ► oo. This discrepancy is due to the finite time necessary to 
establish the constant height boundary condition and can be accounted for by using 
the linear tip solution. Flow for short times t < tn is complicated by the presence 
of inertia and surface curvature in z. However, we show that some of these flows are 
described by similarity solutions. 



Nomenclature 


The dimensional versions of the dependent and independent parameters used in this 
study are denoted by a prime symbol (') unless otherwise specified. The use of the 
subscripted shorthand notation for differentiation should be clear in context. 

Dimensionless numbers and quantities 
a similarity transform parameter 

A cross-sectional flow area 

b similarity transform parameter 

B Hoffman shift parameter 

Bi constant 

Bo Bond number, ApgH 2 /a 

Boh Bond number based on surface curvature, fpgH 2 /a 

Bohl Bond number based on curvature and column length, fpgHC/a 

Ci constant 

Ca Capillary number, pWj a 

Ca t time-dependent Capillary number, p,H/crt 
f surface curvature function 

F generalized similarity function, F(r)) 

T geometric function 

F a cross-flow area function 

F An cross-flow area sub-function 

Fi banded flow resistance function 

h meniscus height from z-axis, h(z, t) 

hi mean curvature 

j aspect ratio function, 1 + Di/D 2 

k geometric function, <£(1 — (j>/2) 

K numerically determined friction factor coefficient [9] 

C tip location, C(t) 

n number of sides on n-sided regular polygon 
Nd index of refraction 

Oh Ohnesorge number, p,j[a pH) l f 2 

P pressure 

P Q scaled pressure, P Q / cos 2 Q 

V wetted perimeter 

q constant volumetric flow rate 

Q volumetric flow rate, Q(z, t) 

7 Z e 2 sin 4 a/ fOh 2 

Re Reynolds number, WH/v 
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5 parameterized surface elevation, S(y , t) 

t time 

T geometric time constant 

u z-component of velocity 

v ^-component of velocity 

w ^-component of velocity 

x spatial coordinate 

y spatial coordinate 

Vm y-coordinate value at contact line 

z spatial coordinate 

Greek Symbols 
a corner half-angle 

(3 numerically determined dimensionless flow resistance [11] 

8 surface curvature angle, tt/2 — a — 6 

e fluid column slenderness ratio 

7 refraction angle 

T slope function for linear tip solution, r(e) 

77 generalized similarity parameter 

rj rr surface viscosity [11] 

6 contact (wetting) angle 

k dimensionless drop volume 

A invariant transform coefficient 

jjl dynamic viscosity 

v kinematic viscosity 

f wavenumber 

p density 

<7 surface tension 

t transformed time (= Fit/ 2) 

<j) 6/a 

Q complimentary corner angle, tt/2 — a 

Dimensional Quantities 
c linear fit coefficient-intercept 

D, Di face width of test cell 
g acceleration field strength, gravity 

gi component of g in i-direction 

g 0 normal gravity, 9.8m/s 2 

G flow coefficient 

iii 



H characteristic meniscus height 

L characteristic length of fluid column, or cell length 

m linear fit coefficient-slope 

R radius of curvature 

R c local mean radius of curvature 

tf curve fit time origin 

t Q inertial reorientation time, response time 1-g to low-g 
W characteristic velocity 

Z global meniscus location on cell centerline 

Vectors 

k inward normal to container wall 
n outward normal to surface S 
ti streamwise unit tangent on S 
t 2 cross-flow unit tangent on S 
v velocity 

Subscripts, Other 

+ invariant transform quantity 

implies dummy variable 
{ ) area-averaged quantity 

0, 1, ... counters for order of expanded quantity, constants, etc. 

I referring to limit a 2 <C 1 

II refers to limit tt 2 <C 1, free surface condition 

II c refers to limit O 2 <C 1, captive surface condition 

cl contact line quantity 

d dynamic quantity 

exp experimentally measured quantity 

H quantity of constant height location 

hyd hydraulic diameter quantity, free surface 
hyd c hydraulic diameter quantity, captive surface 
i counter 

ic inscribed circle quantity 

n-poly quantity for n-sided regular polygonal cross-section 
red quantity for rectangular cross-section 

Tl selected test, 2.0cs PDMS in 22.6mm eq. triangular cell 

T2 selected test, lO.Ocs PDMS in 12.0mm eq. triangular cell 
tip quantity evaluated at tip 

stat static value 


IV 



vise quantity attributed to viscous effect 

wall evaluated at wall-corner axis location of vessel 


v 



1 


Contents 


Preliminary Pages 

Abstract i 

Nomenclature 

List of Tables V1 

List of Figures 

1 Introduction 1 

1.1 Mechanism and Motivation 1 

1.2 Background and Scope of Research 3 

2 Analysis: Governing Equations 6 

2.1 Analysis of Flow in an Isolated Corner 7 

2.2 Asymptotic Equations * • • H 

2.3 Solutions for F* 16 

2.4 Steady Solution 25 

2.5 Infinite Column Solution * • 25 

2.6 Similarity Solutions * • • 27 

2.6.1 Capillary Rise (Constant Height) 28 

2.6.2 Spreading Drop (Constant Volume) 29 

2.6.3 Constant Flow 30 

2.6.4 Generalized Similarity Solution 32 


vi 



2.6.5 Inclusion of Gravity in ^-direction 33 

2.6.6 z-Curvature Consideration 34 

2.6.7 Summary of Closed Form Results (dimensional) 36 

2.6.8 Determination of if, Capillary Rise 38 

2.7 Quantified Constraints 39 

2.7.1 Assumptions for Capillary Rise 40 

2.7.2 Assumptions for Spreading Drop 42 

2.7.3 Assumptions for Constant Flow 44 

2.7.4 Abbreviated Summary of Assumptions 45 

3 Capillary Rise Experiments 47 

3.1 Introduction to Low-g Experimentation 48 

3.2 Description of the Experiment 48 

3.3 Preliminary Measurements 51 

3.3.1 Tip Location, C(t) 52 

3.3.2 Global Meniscus Height, Z(t) 55 

3.4 Detailed Measurements 55 

3.4.1 Surface Profile, h(z,t) for data set T1 57 

3.4.2 Surface Profile, h(z, t) for data set T2 61 

3.5 Discussion 63 

3.5.1 Tip Solution Approximation 68 

3.5.2 Start-up Problem 68 

3.5.3 Further Remarks 71 

4 Conclusion 72 

4.1 Summary of Results 72 

4.1.1 Analytical Contributions 72 

4.1.2 Experimental Contributions 74 

4.2 Recommendations for Future Work 75 

vii 



Acknowledgement 


76 


Bibliography 77 

Appendices 80 

A Analytical Details 80 

A. 1 Equations and Boundary Conditions 80 

A. 2 Fi Details 81 

A. 3 Similarity Details 84 

A. 4 Similarity with g x 85 

B Correction for Optical Distortion 87 


viii 



List of Tables 

2.1 Nondimensionalized dependent and independent variables 8 

2.2 Values of /3, and Fh y d 21 

2.3 Values of K and Fi 22 

2.4 Summarized values from the similarity solutions 31 

2.5 Selected values from the generalized similarity equation 32 

2.6 Summary of constraints 46 

3.1 Test fluids and properties 50 

3.2 Test cell data for capillary rise tests 50 

3.3 Time constants for T1 and T2 for capillary rise 68 

3.4 Time constants for T1 and T2 for start-up problem. 71 



List of Figures 


1.1 Partially filled square container 2 

1.2 Capillary rise due to a step reduction in g 3 

2.1 Sketch of fluid in isolated corner 7 

2.2 Sketch to show variables used in the cross-flow formulation . 13 

2.3 The parameters f and F \ as a function of 8 for various a 14 

2.4 Fi (a, 0) for the free surface condition 23 

2.5 Fi(a,0) for the captive surface condition 24 

2.6 Geometric time constant T(a, 8) 26 

2.7 Similarity function F + (r} + ) for capillary rise, spreading drop, and con- 
stant flow 31 

2.8 Similarity function F + (r] + ) for a > —1/3. • ■ 33 

2.9 Similarity function F(rj) for capillary rise with gravity, g x 35 

3.1 Selected frames from cinefilm records of capillary rise in equilateral 

triangles. 49 

3.2 Cross sections of capillary rise test cells . 51 

3.3 £ vs. t 1 / 2 in acrylic triangle: Effect of viscosity. 53 

3.4 L vs. t 1 / 2 : Effect of cell size and corner angle 53 

3.5 C vs. t 1 ! 2 in acrylic rectangular cells: Effect of aspect ratio. ...... 54 

3.6 Z vs. t 1 / 2 in acrylic triangle: Effect of viscosity. 55 

3.7 Z vs. t 1 / 2 in acrylic cells: Effect of cell shape and size 56 

3.8 C vs. t in acrylic triangular cells with c + m(t — t/) 1 ^ 2 fits 57 


x 



3.9 Z vs. t for T1 and T2 58 

3.10 Traces of the surface height h(z,t) for T1 59 

3.11 h(zi,t) vs. t for T1 60 

3.12 h(zi,t) vs. t for Tl, t > 0.7s 61 

3.13 h vs. z/t 1 / 2 and z/(t — t#-) 1 / 2 for several traces. Data of Tl 62 

3.14 h vs. z/(t — t Q ) x / 2 for Tl with t a = 0.09s 63 

3.15 Traces of surface height h(z,t) for T2 64 

3.16 h(zi,t) vs. t for T2 65 

3.17 h vs. z/t 1 / 2 and zj[t — tn) 1 ^ 2 for several traces. Data of T2 66 

3.18 h vs. z/(t — to) 1 / 2 for T2 with t Q ~ 0.03s 67 

3.19 h/t 1 ^ vs. z/t z ^ for T2 69 

3.20 h(zi, t) as a function of t for T2, Zi ~ zh 70 

B.l Schematic of ray trace correction for h(z,t) measurements 88 




Chapter 1 


We commonly think of capillary phenomena as relevant to small systems such as 
capillary tubes, candle wicks, paper towels, etc. This is because on earth only small 
systems allow surface tension and wetting phenomena to play a dominant role over 
the otherwise controlling influence of gravity. In the presence of gravity the degree 
to which capillary forces can be expected to be significant is measured by the Bond 
number Bo = ApgH 2 /cr, where A p is the density difference across the fluid interface, 
g is the acceleration of gravity, H is a characteristic length scale ( i.e . tube radius, 
pore size, etc.), and a is the interfacial, or surface, tension. The Bond number is 
a ratio of gravitational forces to those of surface tension, and low values, Bo <C 1, 
imply capillary dominated systems. The properties A p and o are essentially fixed 
for a given system. For earth bound systems g is also fixed and capillary dominated 
systems arise only if geometries are small. However, for the low-gravity environment 
of space, surface tension can dominate fluid systems of large extent. 


1.1 Mechanism and Motivation 

The subject of capillary driven flows in containers with interior corners is perhaps 
best introduced via an illustrative example. Figure 1.1 depicts a container with a 
square cross-section partially filled with a liquid in a gravitational field. The majority 
of the interface is flat owing to the large Bond number of the system. However, as 
is commonly observed, especially in the corner regions of the container, the interface 
curves in order to satisfy the contact angle condition along the perimeter of the 
interface. 

A magnified view of the corner region shows that the local radius of curvature 
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Ri of the meniscus decreases as the corner is approached. Since the pressure drop 
across the meniscus is inversely proportional to Ri it follows that the pressure in 
the liquid also decreases as the tip of the meniscus in the corner is approached. 
Thus, a pressure gradient along the corner is established in the wetting liquid due to 
increasing interface curvature. This gradient is in turn balanced by the hydrostatic 
pressure gradient. If gravity were suddenly “turned off,” the interface curvature 
pressure gradient would no longer be balanced and a capillary driven flow along the 
corners would ensue. The extent to which the fluid rises in the corner during the flow 
process is dependent on the contact angle (0) and corner half-angle (a) of the system. 
In many cases the fluid will climb to a great height if allowed, provided that the 
contact angle is low enough. The criterion for when such large scale flows must occur 
is thoroughly addressed by Concus and Finn [l]-[3] who show that capillary driven 
flows in corners of infinite extent proceed to infinite distances when the condition 

6 < tt/2 — a (1.1) 

is satisfied, hereafter referred to as the Concus-Finn condition. 1 This condition is 
satisfied in the case of Figure 1.2 which illustrates this phenomenon for a square 

1 In their mathematical study, Concus and Finn were the first to obtain the result of the discon- 
tinuous behavior of an interface in a corner at the critical value 9 — k/2 — a. The condition of eq 
1.1 appears also in other heuristic studies concerning interfacial phenomena. 
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Figure 1.2: Capillary rise due to a step reduction in g. 

cross-sectioned container employed in a simple drop tower experiment. (Note that 
the figure displays only half of the vessel which is bisected along the diagonal) After 
release of the container in the drop tower, slender columns of fluid rise in the corners 
due to capillary forces and act to drain the container by redistributing the fluid along 
the corners. 

Capillary pumping of this general sort arises in numerous contexts and is particu- 
larly prevalent in low gravity environments where system designs are heavily impacted 
by the presence of container irregularities such as corners. Examples of such systems 
involve most, if not all, in-space fluids management processes, from the positioning, 
control, and transport of liquids, such as fuels in storage tanks, to thermal systems 
such as heat pipes and capillary pumped loops, to the storage and handling of bio- 
logical fluids and wastes. Examples of terrestrial applications include flows in porous 
media and the wetting and spreading of fluids on irregular surfaces-flow processes 
which are commonplace in nature and industry. 


1.2 Background and Scope of Research 

A number of investigations have produced useful results on interface configurations 
(statics), flows, and stability for the capillary behavior of liquids in containers with 
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interior corners. 2 The case of capillary statics has been well addressed by the liter- 
ature. Included in these studies, for example, are the effects of contact angle and 
corner half-angle [1][2], container aspect ratio and corner rounding [4], differing con- 
tact angles on adjacent faces of the corner [5], irregular triangular sections [6], and a 
novel solution approach based on the inclusion of molecular forces [7]. The works of 
Concus and Finn, who have obtained specific results for a host of cylinders of general 
section, in part summarized in [8J, are exemplary and a review of their work may 
serve well the designer requiring insight into low-gravity interface shapes. 

A significant amount of attention has also been devoted to dynamics. Though 
the majority of work [9]- [17] has centered on forced flows and/or wicking (imbibition) 
in “micro-sized” geometries for earth-based applications, low-gravity experiments by 
Masica (discussed in [1] and [8]) and by Langbein [18], in addition to a multitude 
of applied space systems studies [19] [20] have clearly demonstrated the large-scale 
nature of such flows in the absence of gravity. Nonetheless, an ample collection of 
quantitative analytical results on related flows applicable to the unusual conditions of 
low-gravity, is yet to be acquired. Questions of practical importance concerning the 
general response time, flow rate, surface profile, and velocity distribution as functions 
of the container geometry and fluid properties-including the contact angle-remain 
to be fully resolved for the unique range of system parameters accessible with the 
reduction of gravity. 

This work reports an analytical and experimental investigation to answer these 
questions for a limited variety of container types and flow scenarios. The analytic por- 
tion of the work begins with a description of the simplified problem of capillary-driven 
flow along an isolated corner and follows with an asymptotic analysis employing the 
lubrication approximation. The velocity scale used in the nondimensionalization of 
the equations correctly captures the geometric effects of corner angle and contact 
angle and reduces a numerically determined geometric function, used by previous 
investigators, ~ 0(1) to oo, to one which is tightly banded and ~ 0(1). The govern- 
ing partial diff erential equation is first solved by a regular perturbation method for 
the case of the infinite corner. The equation is then transformed yielding similarity 
solutions for the cases of capillary rise, a spreading drop, and a constant flow rate 
condition. The impact of gravity, three dimensional surface curvature, inertia, and 
the moving contact line boundary condition are considered. 

A description of the experimental techniques, selection of test results, and dis- 
cussion thereof follows for the special case of capillary rise. It is shown that the 
analytical results, which often yield closed form expressions, are well-suited for use 
in describing such flows, even at small times. Applications may be made to problems 

interior corners are often referred to as “grooves,” “edges,” or “wedges” in the literature. 
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arising in the low-gravity environment. Only systems which satisfy the Concus-Finn 
condition are considered in detail. For an analysis of the case 0 > tt/2 - a see Lang- 
bein [21] who determines both the shape and stability of fluid columns in corners of 
infinite extent. 
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Chapter 2 

Analysis: Governing Equations 


The particular corner flow problem illustrated in Figure 1.2 was recently studied by 
Dong and Chatzis [17] who employ 0.3-0.5mm diameter tubes of square cross-section 
in an analysis and experiment on imbibition as relates to flows in porous media. As 
in previous and related studies [11] [12], it is sufficient for the authors to simply state 
the assumptions (locally parallel flow, negligible corner axis curvature and inertia, 
etc.) and proceed with solutions to the reduced set of governing equations. The 
‘small’ container size <, O(lmm) of these experiments provides a firm basis for such 
an approach. For a system of equivalent Bond number in a low-gravity environment, 
however, the assumptions come into question. For example, the Ohnesorge number, 




h 

{a pH) 1 / 2 


( 2 . 1 ) 


where p is the dynamic viscosity of the liquid, or rather l/O/i 2 , is the correct choice 
as a measure of the influence of inertia in capillary driven flows; similarly, as is the 
Reynolds number in forced flows. For ground-based studies ( g = 1 g 0 ), a strong 
capillary presence is maintained when the characteristic dimension H scales as the 
capillary length, H ~ (a /pg) 1 / 2 . Therefore, 


1 ap 



1/2 


Thus, in a microgravity environment (g = 10 ~ 6 g 0 ) the parametric range of 1/Oh 2 
is extended up to 3 orders of magnitude. More dramatic yet are the changes which 
might be expected in the viscous time scale t V i SC ~ pH 2 / p. With a 1000-fold increase 
of if, realizable in a low-gravity environment, t visc increases by 10 6 over its normal 
gravity counterpart! To accommodate for such possible scenarios an effort here is 
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Figure 2.1: Sketch of fluid in isolated corner. 

made towards generality. The governing equations are scaled in such a manner that 
the application of asymptotic techniques leads to simply expressed solutions for a 
variety of corner flow problems. A natural means for assessing the validity of the 
solutions arises and the many constraints of the problems are quantified. 


2.1 Analysis of Flow in an Isolated Corner 

Since the corner flow controls much of the flow throughout the container, it is natural 
to analyze the fluid in this region first. An isolated interior corner of angle 2a 
is depicted in Figure 2.1. The corner is partially filled with a fluid satisfying the 
Concus-Finn condition. The coordinate axes are labeled on the figure, the /-axis 
corresponding to the corner axis, and S' is the parameterized height of the free 
surface as measured from the y'-z' plane. 

N ondimensionalization 


As shown in Figure 2.1, L is the characteristic length of the column of fluid along 
the /-axis and H is the characteristic height of the meniscus in the corner along the 
/-axis. The ratio of these length scales (e = H/L) is crucial to the analysis to follow. 
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Lengths 

Velocities 

Other 

* = x’/H 
y = y' /H tana 
z = z'/L 
S = S'/H 
h = h'/H 
£ = CjL 

ym = y' m ax/ Ht ^ a 

u — v! /eW 
v = *// eW tan a 
w = w'/W 
(w) = {wY/W 
W = ea sin 2 a/ (if 

P = HfP’/a 
t = Wt'/L 
A - A! /H 2 tan a 
Q = Q'/WH 2 tana 


Table 2.1: Nondimensionalized dependent and independent variables. 


Table 2.1 provides the quantities used to nondimensionalize the governing equa- 
tions. Primes are used to denote dimensional dependent and independent variables. 
The length scales arise naturally from the geometry of the problem. The pressure is 
scaled using a/Hf, and / = /(a, 0) is a geometric function describing the curvature 
of the meniscus in the x'-y' plane, to be defined later. For this curvature scale to be 
valid the curvature in the x'-z' plane (~ d 2 S'/dz' 2 ) must be small by comparison. 
In other words, the constraint of e 2 / <C 1 must be imposed. The characteristic z'- 
component of velocity, W, is determined through a balance of pressure and viscous 
forces and uniquely incorporates the geometric influence of a on the solution domain. 
The scales for v! and v' are determined via continuity. A passive overlaying fluid is 
assumed and, for the time being, body forces are ignored. 

On substitution of the Table 2.1 quantities, the individual component equations 
of the nondimensionalized Navier-Stokes equation become 


with 


and 


o __ Du 

‘ k dI 

== —P x + e 2 V 2 u 

(2.2) 

9 __ 9 DV 

e 2 1Z tan 2 a— 

i/6 

~ — Py T e 2 tan 2 aV 2 v 

(2.3) 

_ Dw 
11 Dt 

— ~P Z + V 2 w 

(2.4) 

D d d d d 

7 TT = ^7 + u— + v— + w~ 
Dt dt ox ay oz 

(2.5) 

. , d 2 2 d 2 2 . 2 d 2 

= Sim airT + cos^ a . . + e sm a— 
dx 2 ay 2 oz 2 

(2.6) 


« 2 • 4 

(2.7) 

fOh 2 
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Oh is the Ohnesorge number given by Oh = (Ca/Re) 1 / 2 = pf(paH) 1 ^ 2 , where Ca = 
pW/cr and Re = pWH/p are the Capillary and Reynolds numbers, respectively. 
Defined as such, Oh serves as a ratio of the restoring-force time scale of surface 
tension (~ (pH 3 /a) 1 / 2 ) to the viscous time scale (~ H 2 /v ). For Oh 2 <C 1 systems 
are underdamped to ‘abrupt’ disturbances while for Oh 2 1 viscous forces dominate. 
The accessible parametric range of 1Z is increased several orders of magnitude when 
gravity is low. This influence is felt through Oh as described earlier. In general, the 
condition Oh 2 > 1 is not often satisfied even for ‘small’ normal-gravity systems. 

The continuity equation is 

V • v = 0 (2.8) 

The boundary conditions for the above system are as follows: The no-slip condition 
is 

v = 0 (2.9) 

and applies along the walls. The scaled (xH/pW) streamwise and transverse zero 
shear stress conditions on S are 


[(1 - S 2 cot 2 a + e 2 S 2 ) {u y +)v x tan 2 a - S z (w y + e 2 v z tan 2 a) ) 


+2 Sy (U x ~ 

Vy Sz ( 

w x + £ 2 M 2 ) - e 2 Sf ( w z + Uj,))] = 0 

(2.10) 

and 




k x k 2 [(l - e 2 5 2 ) (w. 

x + e? ' u z'j 

- S y (w y cot 2 a + £ 2 v 2 ) + 2 e 2 S z (u x 

~Wz) 

respectively, where 

SSySz 

(uy cot 2 a + v*)] = 0 

(2.11) 

h = 

=(i+^r i/2 

(2.12) 


?r 

tsj 

II 

:(l+|V5| 2 )' 1/2 

(2.13) 

and 

|VS| 2 

— S 2 cot 2 a + e 2 5 2 

y z 



The normal stress condition on S is 


- P + 2 e 2 k 2 2 [M] sin 2 a = 2H (2.14) 


where 

[M] = u x — S y cot 2 a(u y + v x tan 2 a) — S z (w x + e 2 u z ) + S y S z cot 2 a (w y + e 2 v z tan 2 a) 
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(2.15) 


+ SyV y cot 2 a + e 2 S 2 z w z 

and 2 7i is twice the mean curvature of the interface given by 


27 i. = /&2 cot 2 ol ^$yy(l £ 2 S 2 ) + e' 

2 S ZZ ( tm 2 a + S. D - 2 e 2 S y S z S yz 

(2.16) 

Two (of four) boundary conditions on , 
symmetry conditions on y == 0; namely 

S are the meniscus centerline location and 

S = h 

on y — 0 

(2.17) 

f = 0 

oy 

on y = 0 

(2.18) 


where h = h(z,t). The remaining two boundary conditions concern the treatment of 
the moving contact line. This poses particular difficulties in that the physics of the 
moving contact line is not fully understood due to a well known stress singularity at 
the contact line itself [22] [23]. 

The simplest conceivable formulation of the moving contact line boundary condi- 
tion is the case of fluid statics where the condition 

n • k = cos 0 (2.19) 


is applied. 0 is the static (or perhaps equilibrium) contact angle of the fluid-solid 
pair, 1 k is the inward normal to the container walls (on y = ±x) given by k = 
(sin a, =F cos a, 0), and n is the outward unit normal to 5, 

n = (1 + |VS| 2 )- 1/2 (1, -S y cot a, -eS z ) (2.20) 


This boundary condition at the contact line is strictly valid when v = 0 and may 
be approximately correct when the fluid velocity perpendicular to the contact line is 
small. This is, in general, the case for the corner flow problem discussed here where 
the predominant flow direction is parallel to the contact line. The use of eq 2.19 
will be discussed in subsequent sections. The net effect of its application is that the 
two additional boundary conditions for 5, namely, the slope and the location of the 
interface at the contact line, may be determined in terms of h. These conditions are 
given by 

1 + S y cot 2 a _ cos 0 
(1 + |V5| 2 ) 1/2 “ 


sina 


on y = y m 


( 2 . 21 ) 


x If contact angle hysteresis is present and the contact line is in motion, 0 might best be taken 
as the advancing angle on first approximation. For a clear definition of the distinction between 
static and equilibrium contact angle and for a discussion of the difficulties faced when contact angle 
hysteresis is present in systems exhibiting partial wetting see [23, p. 328]. 
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and 


S = y m on y = y m (2.22) 

respectively, the latter condition being the interface /wall intersection condition and 
y m being the value of y at the contact line. 

A mass balance in the ^-direction gives 


BA BQ d i . 

« - - T, ~ 


(2.23) 


where A is the cross-sectional area of the flow in the x-y plane and Q is the volumetric 
flow rate in the ^-direction, namely 


A = 2 f ym (S-y)dy 
Jo 

rym rS 

Q — 2 / wdxdy 

JO Jy 

and { w ) = Q/A. The inlet and outlet conditions at z = 0 and z = 1 and the initial 
conditions remain to be addressed. 


2.2 Asymptotic Equations 

Though eqs 2, 2-2. 4 and their associated boundary conditions eqs 2.9-2.22 are a 
formidable set of equations, the slenderness ratio e, appears as an obvious choice 
for a small parameter upon which to base an asymptotic analysis. The limit e 2 <l 
is the common lubrication approximation and, as in many other contexts [24], is 
motivated by observations of the flow phenomena such as that depicted in Figure 
1 . 2 . 

For the case e 2 1 the dependent variables of the problem may be expanded in 
asymptotic series 


u = u 0 + e 2 Ui + ... 

(2.24) 

v = v 0 + e 2 Vi -f ... 

(2.25) 

w = w 0 4- e 2 wi + ... 

(2.26) 

P = P Q + e 2 P\ + ... 

(2.27) 

S ^ So + € 2 5 1 + ... 

(2.28) 
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where terms of 0(e) are absorbed into the leading order terms for convenience. In 
the limit e 2 <C 1, assuming for the present 7Z ~ 0(e 2 ), the leading order equations 
reduce to the ^-component momentum equation 


dPo 

dz 


d 2 w 0 2 d 2 w„ 
= sin* ol ^ n + cos a- 


dx 2 


dy 2 


(2.29) 


Also to leading order, the boundary conditions for eq 2.29, eqs 2.9, 2.11, and the 
constraint of flow symmetry about y — 0 may be expressed as 


w o = 0 on y ~ x 


dw Q dS 0 dw Q 


dx 


dy dy 
dw 0 
dy 


cot 2 a = 0 on x — S Q 


= 0 on y = 0 


respectively. The normal stress condition can now be written as 

2 \ -V2 


P 0 = /cot 2 a 


d 2 S 0 

dy 2 


1 + 


'dSoY 

, d v ) 


cot 2 a 


with the leading order boundary conditions 

S 0 = h on y = 0 


dy 
555 ym 

1 j~ ^dy C°t 2 a 


0 on y = 0 


( i + c ° t2 a ) 


\ 1/2 


on y = y m 

cos 0 
sin a 


on y = y m 


(2.30) 


Noting that P 0 — P 0 (z,t ), the normal-stress equation may be solved directly with 
its associated boundary conditions to find S Q , y m , /, and P 0 (h). This procedure 
results in forms of these quantities which are unnecessarily awkward. Less cumber- 
some forms are obtained by solving the geometrical ‘cross-flow’ problem sketched in 
Figure 2.2 for the pressure in the liquid as a function of meniscus height h, with the 
dimensional radius of curvature 2 taken to be constant R' = fh f . This yields 

Po = -\ (2.31) 

2 Note that the curvature of the interface is not a constant in the dimensionless x-y plane. 
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ym 


Figure 2.2: Sketch to show variables used in the cross-flow formulation, eqs 231-234. 
The height of the interface at a given z is S 0 (K y), where h{z,t) is the height at 

y- 0 . 


s 0 = Mi + f)~f h 



2/ tan as \ 

fh ) 


1/2 


fh sin 8 
tana 


/ cos 0 \ ' 1 / sin(a + 6) 

\sina ) \ sin a 



(2.32) 

(2.33) 

(2.34) 


where 8 = 7r/2 — a - 6. The parameter <5 is introduced as a preferred angle (as 
opposed to 6) which measures the degree of curvature of the interface in the x-y 
plane. The condition 0 < 8 < tt/ 2 is the Concus-Finn condition. 3 The use of / in the 
nondimensionalization of the problem leads to the unique expression for P 0 given in 
eq 2.31 and the sign of f indicates the anticipated flow direction of the z-component 
velocity scale W. The parameter / characterizes the strength of the driving force for 
the flow due to curvature of the interface, independent of the cross-flow area. The 
parameter Fa = F A (a , 8) is the cross-sectional area function {A = h ,2 F A ) and may 


3 The case of 8 < 0 is also accounted for in this result. Note that / can change sign depending 
on the value of <5 indicating a change in positive to negative curvature of the interface. Positive 
curvature is depicted in Fig. 2.2. 
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Figure 2.3: The parameters / and Fa as a function of 8 for various a (deg). 


be determined geometrically 


F a 


f 2 sm 2 8 ( (2 8 — sin 28 \ 

1 — — rr tan a 

tan a \ \ 1 — cos 28 ) 


(2.35) 


Figure 2.3 shows Fa and / plotted against 8 for a variety of corner half-angles a. 
Note that Fa is an 0(1) quantity for intermediate values of a (tt/12 <, a < 57r/12) 
and shows weak functional dependence on 8 for fixed a. It can be shown further 
that 1 < F^/tana 4/3 for all a and 8 , and thus, the length scale y' ~ H'tana 
is appropriate for all a and 8 and that, in general, Fa ~ tana. Note also that / 
behaves like 1/8 for <5 — > 0. 

Regardless of the technique employed to solve eq 2.29 the resulting form for (w Q ) 
will be 

{Wo) = - Fih >f = - Fi f z (2.36) 

where Fi is some positive geometric function of a and 8 . It is interesting that viscous 
resistance (~ hr 2 ) is balanced by capillary forces (~ h~ 2 h z ) in such a way as to leave 
(■ w 0 ) only dependent on the slope of the interface h Z) and not on h. As a result, 
non-zero velocities are possible for h = 0 provided h z ^ 0. 
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Noting that A oc h 2 , eq 2.36 may be substituted into eq 2.23 which, after some re- 
arrangement, yields the leading order governing equation for the corner flow problem 
in the e 2 <l and fe 2 <C 1 limits 





(2.37) 


In this way, the meniscus centerline height h becomes the principal dependent variable 
of the problem. Eq 2.37 also appears in the literature where it describes unsteady, 
nonlinear heat conduction of the form T t = (T 1/2 T 2 )^ [17] [25]. 

In addition to the selection of boundary conditions, solutions to eq 2.37 need only 
Fi to begin. Generally speaking, Ransohoff and Radke [11] and Ayyaswamy et al 
[9] have determined Fi for this problem. 4 They first use a finite element scheme to 
solve eq 2.29 for w Q . They then integrate over the cross sectional area to determine 
(w 0 ) and the tabulated results are presented in terms of either a dimensionless flow 
resistance, /?, in the case of [11], or a dimensionless friction factor coefficient, K, as 
in the case of [9]. Ransohof and Radke [12] also employ eq 2.37, as do Dong and 
Chatzis [17]. Both employ fi in its solution. In [12] the continuity equation in the 
general form of eq 2.37 is solved numerically for a problem relating to snap-off in 
constricted pores while in [17] a similarity solution is found applicable to imbibition 
in square capillary tubes. The usefulness of these solutions is that they are only 
restricted to e 2 < 1 (slender fluid column), e 2 f < 1 (slight curvature along 2 :- axis), 
and 1Z <, 0(e 2 ) (small inertia). For many practical systems / ^ 0(1) and the second 
constraint is naturally satisfied by the first. 

An analytic form for Fi, however, is still desirable and under certain restricted 
cases further asymptotic analyses are possible to determine its functional form and 
hence the leading order behavior of the corner flow, namely h(z , t). It will be demon- 
strated that the scaling presented here renders Fi a weak function of a and 8 such 
that 1/8 <, Fi < 1/6 for all values of a and 8 . The utility of this fact should become 
apparent in subsequent sections where problems of practical interest are addressed. 


4 Lenormand and Zarcone [26] and Trefethen [27] use a hydraulic diameter approach in effect to 
determine Fi, the utility of which is discussed in [11]. Lenormand and Zarcone ([26], Appendix 2) 
also anticipate the work of [17]. 
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2.3 Solutions for Fi 


In this section asymptotic techniques axe applied to eq 2.29 for the limiting cases 
of small (I) and large (II) corner angle. 5 The latter is solved for the free (II) and 
captive (II C ) surface conditions. A form for Fi is also determined using the hydraulic 
diameter approach (III, F hyd ). In addition, two exact analytic expressions for 
are possible for the special cases of 0 = a = tt/ 4 and for 9 = 7t/3, a = 7r/6. These 
results are then compared to the numerical results of Ransohoff and Radke [11] who 
determine (5 for a wide range of a, 0 , and surface conditions satisfying the Concus- 
Finn condition. From the definition of /?, it can be shown that 


Fi 


/ 2 

p sin 2 a 


(2.38) 


For the free surface problem, the numerical results of Ayyaswamy et al [9], who 
calculate a friction factor coefficient K for the problem of gravity driven flow in 
triangular grooves, are also applicable and are used to compare with F{. It can be 
shown that 

- (— ) 

K \f sin.6 J 


F< = 


(2.39) 


An outline of the general asymptotic solution is provided for case I only. The 
analytical details of the solutions for Fu, F Ud and Fhyd are relegated to Appendix 
A.2. Though only the zeroeth order terms are applicable due to the nature of eq 
2.29, higher order terms are retained for completeness. 


Case I . Fj, Small corner angle solution, tan 2 a < 1 
In the limit a 2 <C 1, an expansion for w Q may be written 

w 0 = w 0o + w Ql tan 2 a + w 02 tan 4 a + ... (2.40) 


and when substituted into eq 2.29 may be solved for w 0 to any desired order subject 
to the no-slip and symmetry conditions. This procedure yields 


1 


W n = 


2 cos 2a dz 


&F 0/2 2 \ 

(y - * *) 


(2.41) 


The average velocity, ( w Q )i , may in turn be determined using 

tana dP 0 rSo 


( ™ 0 ) i = 


h 2 F A cos2a dz 


n s 0 

(y - r/)dxdy 


(2.42) 


5 These cases were anticipated in [11] yet no analytical treatment of them has appeared in the 
literature to my knowledge. 
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Upon substitution of the quantities P 0 , 5 0 , y m and / from eqs 2.31-2.34 the above 
integration yields 

K)/ = -f/f; (2-43) 

where Fj replaces Fi of eq 2.37 for the limiting case of small corner angle, tan 2 a<l, 
and is given by 


Fj = 


2/ [sin <5 / 


Fa cos 2a [6 4 

sin 28' 


8 — 2 sin 8 + 


sin 2(5 


£ 

2 


V r sin 28 sin 3 <5 
<5-2sm<5-f h 




_5_ 

16 


(5 + 


\ cot 2 o\ 2 . r sin 8 cos 3 8 

'--r - 3 m<+ 


24 


sin 3 5 sin 4 (5 cot 3 a 


6 sin 2 a 


12 


2 ‘ 3 sin 2 aj 

^1 -j~ 3 cot 2 + 

(2.44) 


Some limiting values of Fj are 


^« 1 )-i + f + S +0(dt ) 


Fi(6 = 7r/2 - a) = 1 + ^ 




57T 

12 



+ 0(a 3 ) 


which reveal that in the limit tan 2 a — » 0, Fj — > 1/6, and thus (w 0 ) = — h 2 /6. Note 
also that <5 = n/2 - a when 0 = 0. 


Case II . F/j, Large corner angle solution, a 2 < i, 8 2 < i 

An asymptotic solution for Fi is also possible when a approaches 7r/2. Introducing 
Q, = 7r/2—a, the case of the large corner angle limit is given by Q, 2 1. For simplicity 

of presentation, <j> = 8 /Cl. The limit O 2 <C 1 automatically implies < 1 when the 
Concus-Finn condition is satisfied. 

A. Free Surface Condition. A similar analysis as the above yields the full expression 
for Fi = Fj / to 0(Q 2 ) 


x 



70 — 126<f> + 77<fi 2 — 160 3 , n 2 

11 ~ 35 (</> — 2) 2 (3 — 2<t>) + (0 - 2) 2 (3 - 2<t>f 

4580 449<£ 2 2840 3 16130 4 57^> 5 160 6 \ 4 

+ 21 ^ 700 + "350" ~ "225" J + *■ ) 


(2.45) 
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with the limiting behavior 


i n 2 

FtM = 0) = - + T + o(n 4 ) 

1 4Q 

fV/W-V-^ + j^S! 2 + <>(«•) 

B. Captive Surface Condition. This solution follows similarly for the case fl 2 1, 
but with a captive (no-slip) surface condition. The boundary condition on x = S Q 
becomes w = 0. The full expression for Fi = Fji c to 0(Q 2 ) is 

_ 70 - 1260 + 770 2 - 1 6<£ 3 

^/j c = — — — rr: 1“ 


Q 2 


(<t> — 2) 2 (3 — 2<^>) 2 
with limiting behavior 


100 73 (t> 2 

3 + 30 


140 (0 — 2) 2 (3 — 20) 
340 s 


67 0 4 2 90 5 


40 s 


105 175 175 225 


+0(fl 4 ) (2.46) 


*».(*-«)- 1 + ln 1 + 0(f!‘) 


ft.W'-D-l + iff + oifi*) 

F/7 (0 = -l) = — + 110 - 3 n a + o(n 4 ) 
' 6300 354375 v ' 


Case III . Fhyd , Hydraulic diameter solution 

A. Pree surface condition. Application of the hydraulic diameter concept for this 
problem gives the result 


Fhyd = 


1 

8 


1 -( i 7^) tana 


(2.47) 


with limiting cases 

Fhyd{6 = 0 ) = - 

1 1 ^O 2 

F^«1,0=1) = - + — + 0(0 4 ) 
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F hyd (rf « i, * = -l) = | + ^ + o(n 4 ) 


B. Captive surface condition. For the case of the captive surface condition the hy- 
draulic approach yields the full expression 


Fhyde ~ ~ 


1 

'l- 

( 2<5-sin2<$ 
( 1— cos2£ 

) tana 

2 

1 

8 

1 - 

( 1-cos (5 \ 
V sin 8 ) 

tana 


1 i 8 sin a 
L 1 ~ r sin <5 J 


(2.48) 


with the limiting behavior 


F'hydcifi — 0 ) — 


1 

8(1 + sin a) 2 


1 4QQ 2 

F W ^<1, # = 1)--L + 3 S5 + 0(!1‘) 

O 2 

Fhydcffi C = -1) “ + — + 0(Q 4 ) 


Case IV . Fjv, Exact solution for a free surface, a = 0 — n/A 

The analytic solutions for F { are adapted from classic solutions of viscous flows in 
noncircular ducts [28]. For the case a = 7r/4, 9 = 7t/ 4, <5 = 0, and the free surface 
(zero shear stress) boundary condition it can be shown that 


Fjy = 


1 

3 


1 


192 A tanhz7r/2 
75 


7 r 


3 , 5 ... 


= 0.14057... 


(2.49) 


The actual flow domain here is the half section of a square duct bisected along the 
diagonal. 


Case V . Fy , Exact solution for a captive surface, a = 7r/6, 0 = 7r/3 

This is the solution for viscous flow through a duct with equilateral triangular cross 
section. The boundary condition on the fluid interface is the captive (no slip) surface 
condition [28]. Again, for this case, 5 = 0 and 

F v = (2.50) 

15 


Summary of solutions for F* 
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Table 2.2 summarizes the above results for F* and includes the numerical values 
of [11]. Decimal values for F* listed in the table are values determined using the 
numerical value of j3 and eq 2.38. Fractions or values are either asymptotic or 
exact analytical results. 

Observation of the tabulated quantities reveals that, though (3 yields values from 
0(1) to oo for the full range of 6 and a satisfying the Concus-Finn condition, F* is 
banded and 0(1) particularly for the more important free surface condition where 
1/8 <, F < 1/6. This is confirmation that the governing equations are scaled appro- 
priately by the quantities of Table 2.1 and that sin 2 a and / appearing in W capture 
the respective geometric components of viscous resistance and capillary driving force 
for the flow. For the case of the Concus-Finn condition with a free surface (3 > 6. If 
0 — yr/2 — a, then (3 = oo. 

In all cases Fh y d < Fi and in the extreme cases, Fi can be as much as 40 times 
larger than F hyd (a = tt/2 and 0 = 0), or differ by as little as 10% (a = 0 = 7r/4). 
The hydraulic diameter approach does not nearly capture the nature of the correct 
Fi over the full range of a and 0, and though the prospect of an analytic form for F 
for all a and 0 is attractive, the results/warnings given by Ransohof and Radke [11] 
should be heeded and the hydraulic diameter approach should not be followed. 

Table 2.3 compares Fi to K determined numerically by Ayyaswamy et al. [9] for 
a selection of a and 0. It can be shown from the asymptotic analysis that for a 2 <C 1, 
K — 48 4- 0(a 2 ). In the large angle limit fF « 1, K — (8/F)( 1 - <5/30 + 0(Q8)). 
Note also that if <5 = 0, F = 8/K. In contrast to (3, K(a, 6) is a banded function 
similar to F, and, as seen from Table 2.3, takes values in the range 30 £ K < 57 for 
the values of a and 6 listed. 6 The primary differences between K and F are sourced 
in the length scale choice for the cross-flow problem. Ayyaswamy et al. [9] use the 
distance of the contact line from the corner to scale all lengths whereas h (or H) is 
used in this analysis. The fact that Fi is more tightly banded than is K for the full 
range of a and <5, and that F ~ 0(1) while K ~ 0(10), indicates that H is a more 
correct choice for the cross-sectional length scale of the problem. 

Figure 2.4 displays the results of Tables 2.2 and 2.3 for systems with a free surface. 
It is readily observed that F falls within a narrow band of values for all a and <5. 
In Figure 2.5, F is plotted against a for the captive surface condition. The banded 
behavior of F for this problem is reduced. The asymptotic solutions for (] 2 <1 are 


6 K may yield values less than 30 for a > 60°. For example, using the value of Fi determined 
from the value of /3 for a = 72°, 6 = 0, eq 2.39 may be employed to show K = 27.9 for this 
case. However, F* values determined from the numerical values of (3 and K show a discrepancy of 
approximately 3% between the two numerical results at low contact angle (compare F in Tables 
2.2 and 2.3 for a = 30°). 


20 





Free Surface (t)rr = 0) 

Captive Surface (t]rr = 10 s ) | 

a 

e 

P 

Ft 

Fhyd, 

P 

Fi 

F } iyd 

0 

0 

6* 

1/6 

1/8 

6* 

1/6 

1/8 

10 

0 

9.981 

0.146 

0.0548 

13.6 

0.108 

0.0353 

15 

0 

12.95 

0.141 

0,0382 

19.91 

0.0914 

0.0209 

30 

0 

31.07 

0.129 

0.0154 

65.02 

0.0615 

0.00600 

36 

0 

46.67 

0.127 

0.0115 

108.1 

0.0544 

0.00405 

45 

0 

93.93 

0.124 

0.00790 

248.8 

0.0469 

0.00248 

60 

0 

443.0 

0,126 

0.00498 

1411 

0.0395 

0.00137 

72 

0 

3185 

0.131 

0.00395 

11390 

0.0367 

0.00102 

90 

0 

00 * 

1/7 

1/8 

oo* 

1/28 

1/32 

30 

0 

31.07 

0.129 

0.154 

65.02 

0.0615 

0.00600 

30 

5 

30.89 

0.132 

0.0205 

65.64 

0.0619 

0.00813 

30 

10 

31.94 

0.133 

0.0261 

68.24 

0.0624 

0.0106 

30 

20 

38.10 

0.136 

0.0396 

81.58 

0,0634 

0.0166 

30 

30 

54.09 

0.138 

0.0558 

115.9 

0.0644 

0.0241 

30 

36 

75.20 

0.139 

0.0671 

161.3 

0.0649 

0.0292 

30 

45 

165.2 

0.141 

0.0861 

355.5 

0.0656 

0.0380 

30 

60 

00 * 

- 

1/8 

00 * 

1/15 

1/18 

45 

0 

93.9 

0.124 

0.00790 

248,8 

0.0469 

0.00248 

45 

5 

93.99 

0.127 

0.0127 

253.2 

0.0473 

0.00408 

45 

10 

100.2 

0.129 

0.0191 

270.8 

0.0478 

0.00622 

45 

20 

139.0 

0.133 

0.0373 

374.9 

0.0493 

0.0124 

45 

30 

290.7 

0.136 

0.0640 

782. 

0.0506 

0.0217 

45 

36 

698.2 

0.138 

0.0850 

1881 

0.0512 

0.0291 

45 

45 

oo* 

0.14057... 

1/8 

00 * 

- 

1/(12 + 8\/2) 


Table 2.2: Values of /?, F;, and Fhyd • An asterisk * indicates that the value was 
determined by eq 2.38 using F*. The parameter rj RR is the surface viscosity coefficient 
[11]. The free surface condition uses r) RR = 0, and for these numerical solutions, the 
captive surface, rj RR = oo, uses the approximation rj RR = 10 8 . 
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K 

Fi 

K Fi 

K 

Fi 

e 

a = 

5° 

a = 

10 

a — 


0.1 

45.691 

0.159 

43.784 

0,152 

40.588 


5 

46.295 

0.159 

44.951 

0.152 

42.814 

1 

10 

46.798 

0.159 

45.927 

0.152 

44,686 

0.141 

20 

47.539 

0.159 

47.376 


47.483 

0.142 

30 

48.031 

0.160 

48.354 

0.154 


0.144 

40 

48.356 

0.160 


0.155 


0.145 

50 

48.559 

0.161 

49.454 

0.156 

51.712 

0.147 

60 

48.665 

0.162 

49.714 

0.157 

52.367 

0.149 

70 

48.686 

0.163 

49.824 

0.158 

52.778 

0.151 

80 

48.625 

0.164 

49.793 


- 

- 

85 

48.563 

0.164 

- 

- 

- 

- 


a — 

30 

a = 

40 

a = 

50 

0.1 

37.797 

0,132 

35.215 

0.128 

32.797 

0.126 

5 

41.053 

0.133 

39.544 

0.129 

38.346 

0.128 

10 

43.790 

0.134 

43.149 

0.130 

42.863 

0.130 

20 

47.864 

0.136 

48.432 

0.133 

49.262 

0.134 

30 

50.681 

0.138 

52.058 

0.136 

53.563 

0.138 

40 

52.706 

0.140 

54.703 

0.139 

56.720 

0.141 

50 

54.201 

0.142 

56.762 

0.141 

- 

- 

60 

55.316 

0.144 

- 

- 

- 

- 


a == 

60 



a = 90 

(« = 0) 

0.1 

30.449 

0.127 



48* 

1/6 

5 

37.615 

0.130 



a = 90 

(« = 0) 

10 

43.214 

0.133 



37.3* 

1/7 

20 

50.571 

0.140 



\ a = 0 (for all 9) 

30 

55.315 

0.144 



48* 

1/6 


Table 2.3: Values of K and F* for the free surface problem. An asterisk * indicates 
that the value of K was determined by eq 2.39 using the asymptotically determined 
value of Fi. 
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Figure 2.4: Fi(a,6). For the free surface condition the value of Fi falls between 
the curves for 8 = 0 and 6 = 0 approximated by dashed lines. For all a and <5, 
1/8 <, Fi < 1/6. Data from Tables 2.2 and 2.3. O A O K, A asymptotic, □ Case 
IV. 


extended across the full domain of a. From this point on only systems satisfying the 
free surface condition as well as the Concus-Finn condition will be discussed. 

With Fi determined analytically for cases I or II above, or by noting that 1/8 <, jFJ < 
1/6 is a weak function of a and 0-the precise value of which may be determined nu- 
merically [9][ll]-eq 2.37 can now be treated using asymptotic and/or other analytical 
techniques. Introducing r = Fit/ 2 transforms eq 2.37 to 


dh _ 0 f dh \ 2 ,d 2 h 

dr \ dz ) dz 2 


( 2 . 51 ) 
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Figure 2.5: Fi(a,9) for the captive surface condition. Dashed lines are asymptotic 
results for Q 2 <C 1 for 8 = 0 and 0 = 0. Data from Table 2.2. O P (solid line 0 = 0), 
A asymptotic, □ Case V. 
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2.4 Steady Solution 


The steady solution to eq 2.51 is 

h = (B 1 + B 2 z) l l* (2.52) 

where Bi and B 2 are given by the boundary conditions. By specifying h(zi,r) and 
h( z 2 ,r), Bi and B 2 are determined, and eq 2.36 may be solved and substituted into 
Q = A{w) which gives the maximum flow rate due to capillary pumping for a fluid 
column of length C. 


2.5 Infinite Column Solution 


The principle dependent variable h may be expanded as 

h = h Q + chi + c 2 h 2 T (2.53) 


Eq 2.51 includes terms through 0(e). On substitution of eq 2.53 into eq 2.51, standard 
asymptotics yields the 0(1) equation 


dh Q _ ( dh 0 \ 2 d 2 h 0 
dr ~ 2 {dz ) °dz 2 


(2.54) 


and the 0(e) equation 


dhi dh Q dhi d 2 hi d 2 h 0 

-fr = ^-d7 + h °-w +hl -d* 


(2.55) 


One (of many) simple solution may be obtained for the corner flow by noting that 
h Q = const is a solution to eq 2.54. This solution describes an interface at steady 
state in a corner of infinite extent. Then, by solving eq 2.55, h may be determined 
to 0(e) as 

h = h 0 + e exp[—^ 2 h 0 r]Ci cos (£z + C 2 ) + 0(e 2 ) (2.56) 

where $ is the dimensionless wave number (f = £'/ L) of some disturbance and C\ and 
C 2 are constants given by the boundary conditions. Prom this result several useful 
characteristics for the flow may be distinguished. 

A key feature of eq 2.56 is the emergence of a time constant. When redimension- 
alized , the decaying exponent reveals 


£ 2 h 0 r ~ e 


, Fi sin 2 a at 

7 


= e\TCa t ) 


-l 


(2.57) 
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Figure 2.6: Geometric time constant T(a , <5). 

where H is the unperturbed height of the meniscus. T = f /Fi sin 2 a and Ca t — 
Hfi/crt , the latter being a time-dependent Capillary number. 

Recalling that Fi & 1/7, for given base state fluid height H and fluid properties 
the decay function will be a t sin 2 a/f. Thus, the time constant will be oc 
//sin 2 a. The relationship T = //sin 2 a is plotted in Figure 2.6 for a variety of 
a and 8 , for cases satisfying the Concus-Finn condition. As can be seen from the 
figure, the time constant rapidly increases as <5 — > 0, particularly for 8 ^ 10° where 
T ~ 1/(8 sin 2a). The lowest value for T is achieved for a = 30°, 8 = it/ 2 — a. This 
remains true when the slight variation of Fi with a and 8 is considered. The function 
T is the ratio of geometric quantities which characterize the viscous (~ 1/ sin 2 a) and 
capillary forces (~ 1//). Changes in fluid behavior due to contact angle are effected 
solely through the latter. 

The contact line velocity can be determined from eq 2.56. The associated Cap- 
illary number may also be estimated which shows an initial maximum value of 
Ca c i ~ e 3 /T that decreases exponentially with time. This value of Ca is indeed 
small and supports the use of the static contact angle slip condition at the contact 
line, eq 2.19. 

It can be shown that the inclusion of gravity along the rr-direction ( g x ) modifies 
the result of eq 2.56 to produce 

h = h 0 -F eexp [-£ 2 h 0 (l + h 2 0 Bo H ) t] C\ cos (£z + C 2 ) (2.58) 
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The form of eq 2.56 is recovered for Bo u = fp9xH 2 /a = 0 (see Appendix A.4) and 
the column is unstable for Bo H < -1 (h 0 = 1). Order of magnitude predictions of 
settling times for slightly disturbed systems are possible by suitable selection of H 
and L. 


2.6 Similarity Solutions 

Using similarity transformations eq 2.51 may be solved for several problems of prac- 
tical interest. In this section a generalized similarity ODE for h(z , t) is derived. Then 
solutions are found for the special cases of capillary rise, a spreading drop, and con- 
stant flow. All of these cases possess an advancing meniscus where h = 0 at the 
tip. The location of the tip on the z-axis is denoted by C = £(r). These three 
problems were alluded to in [17] and [26] and the problem of capillary rise was solved 
in the former. The special case of capillary rise is rederived along with the other new 
solutions. 


The introduction of 

h = C X T a F{ rj) t) = C 2 zr b C = riti P C 2 1 r- b (2.59) 


where Ci, C 2 , and r) tip are constants, transforms the governing eq 2.51 to 

_l-a-26 

FF nn + 2F] - [aF + br,F v ] = 0 (2.60) 

The case of a = b = 0 is the steady state problem which yields a simple analytical 
solution, eq 2.52. Similarity in eq 2.60 is achieved when 


6 = _( 1 ± 2 ) 

2 


(2.61) 


giving 

FF m + 2 F„ 2 - ^ [aF + br,F n \ = 0 (2.62) 

The difficulty now faced is one of identifying practical boundary conditions which 
yield to the similarity form and supply values for a and b. The constants C\ and C 2 
may be selected as desired to simplify eq 2.62. For the special cases to be introduced 
below, the derivation and/or selection of the quantities a, 6, C\ and C 2 is shown in 
Appendix A.3 (see also [25]). 
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It is useful to note that, in the case of similarity, eq 2.62 is invariant under the 
transformation F = X 2 F + and rj = A rj + . The common requirement of F(rj tip ) = 
F + {r}tip) — 0 for problems involving an advancing tip permits the determination of 
Fj+inip) from eq 2.62, namely 


= (2 - 63) 

Due to the added degree of freedom introduced by A, may be conveniently set 
equal to 1. Eq 2.63 then specifies the slope at the tip. Now, knowing both F + and 
F++ at the tip permits the use of a single step backwards Runge-Kutta method, when 
numerical solution is necessary. The parameter A, and hence F and 77, may then be 
determined by the integral mass balance or the applicable boundary condition. 

It is remarkable to observe from these equations that the slope condition at the 
tip is given by the equation itself. Though the magnitudes of the nonlinear terms 
2 F++ 2 and rj + F++ are maximum at the tip, the linear terms contain F + which is 
approaching zero there. As a consequence, F p+ divides out and the result is a linear 
system near and at the tip requiring only satisfaction of F + (l) = 0. 

The invariant similarity equations for the three flow scenarios are presented below 
with their respective boundary conditions choosing rj^ p = 1. 

2.6.1 Capillary Rise (Constant Height) 

This is the problem of an initially flat interface (Bo ~ 00) in the x~y plane. Gravity, 
acting in the negative ^-direction, is suddenly “turned off 5 and the flow proceeds up 
the corner in the positive ^-direction due to capillary forces-a process comparable 
to the phenomena depicted in Figure 1.2. For this problem a = 0, b = —1/2, and 
C\ = 1, C2 — 2 -1 / 2 are selected. The choice of a and b result from an assumption of 
constant meniscus height at z = 0, or rather, h(0, r) = 1, which leads to F(0) = 1 
for the transformed problem. 

The condition F(0) — 1 is not obvious as the appropriate choice for this ‘capillary 
rise’ problem. Though a condition of this type has been used previously [17], it has 
not been demonstrated that in a time-dependent flow the condition .F(O) = 1 is 
established at rj = 0. The validity of this boundary condition will be discussed in 
conjunction with the presentation of experimental data acquired to this purpose. 
An analytic means for determining the constant height for several simple cell cross- 
sections is provided in a following Section 2.6.8. [Note that h = F, 77 = ^(2r)" 1 / 2 , 
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(2.64) 


C = T] t ip(2r ) 1 / 2 .] The invariant eq 2.62 becomes 

F + F* +v+ + 2 Ftf + 7 fF$ = 0 

subject to 

F + ( 1) = 0 Ffc( 1) = ~ 

Numerical solution to the above system gives F + (0) = 0.345... from which A may be 
determined in order to satisfy F(0) = 1. Quantities of practical interest, namely, the 
location and slope of the interface at the tip, and the slope of the interface at r\ = 0 
may be determined as 

V tip = A = (F + (0))' 1/2 = 1.702... (2.65) 

Fr,(ritip) = -A/2 = -0.851... 

F n ( 0) = AF+(0) = -0.349... 

The numerical solution for F + (rj + ) is shown in Figure 2.7. 

A first order asymptotic solution near the tip, where F + — ► 0 and r} + — * 1, gives 

F + = i (l - J? + ) (2-66) 

and so r) tip = 2 1/2 and ^(0) = F v (r) tip ) = ~2 _1/2 . It is interesting to note that this 
solution, though applicable to the tip region only, also gives F( 0) = h(0) = 1, that 
it satisfies F+ +r}+ = 0 for all 77 + , and that rfap is reduced by approximately 17% when 
compared to rj tip determined by the full solution of eq 2.64. The tip solution of eq 
2.66 is sketched in Figure 2.7. 

2.6.2 Spreading Drop (Constant Volume) 

This is the case for a drop of volume H z being instantaneously placed in an interior 
corner. The drop spreads in both the positive and negative ^-directions. Symmetry 
is maintained across the x-y plane at 2 = 0 and the height of the drop at z — Q 
decreases as the fluid is pumped along the corner. A conservation of mass integral 
gives a = -1/5, b = — 2/5, and C\ = (rc 2 /5) 1/5 , C 2 = (5 2 /c)“ 1/5 are selected, where 
k = z/2Fa- [Note that for this case h = k; 2 / 5 (5t)“ 1 / 5 F, 77 = kT 1 / 5 2:(5t) -2 / 5 , and 
C = « 1/5 ?7tip(5r) 2/5 .] 

The governing invariant eq 2.62 for this problem is 

F + F+ „ + + 2 F+ 2 + 2 V + F+ + F+ = 0 (2.67) 
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subject to 


F + (1) = 0 F+( 1) = -1 

for which an analytic solution exists 

F + = l (i - V +2 ) 


( 2 . 68 ) 


A may be determined via the integral mass balance and quantities of interest may be 
computed, 

— awr-( if—. 


(2.69) 


Fr,(Vti P ) = - (y) = -1.496... 


Eq 2.68 represented in terms of h(z,r) is 

h = 1a 2 ^ 5 (5 t)- 1 ' 5 (l - rf/ A 2 ) 

where 0 < rj < A. The complete result for F + (rj + ) is provided on Figure 2.7. 


(2.70) 


2.6.3 Constant Flow 

When the flow rate q at the inlet (z — 0) is held fixed, the amount of fluid in the 
corner, initially zero, increases linearly with time and the fluid profile increases both 
in height and length. The integral mass balance for this problem gives a = 1/5, 
b = -3/5. Constants C Y = (5<? 2 ) 1/5 and C 2 = (5 3 g)“ 1/5 are chosen. [Note that h = 
g 2/5 (5r) 1/5 F, rj = g~ 1/5 z(5r)“ 3/5 , and C = q 1/5 r] t ip {^ r ) Z/5 , where q = q , /F i WH 2 F A ] 


Eq 2.62 can now be written as 

F+F+ „ + + 2 F+ 2 + 3 - F+ = 0 (2.71) 


subject to 

F+(1) = 0 J$(l) = -§ 

Numerical solution of the above gives F + (0) — 1.200... and F^(0) = 0.935... A may 
be determined using the mass balance integral 


A — Tjtip — 



0.942... 


(2.72) 
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Figure 2.7: Similarity function F + (rj + ) fox the (a) capillary rise, (b) spreading drop, 
and (c) constant flow problems. The constant slope line (d) is the approximation valid 
near the tip for the capillary rise problem. Note that for this solution F + (0) = 1/2. 


Case 

Vtip 

m 

F*(0) 

F v (Vu P ) 

a. Capillary Rise 

1.702... 

1 

-0.349... 

-0.851... 

b. Spreading Drop 

1.496... 

1.119... 

0 

-1.496... 

c. Constant Flow 

0.942... 

1.065... 

-0.881... 

-1.413... 

d. Linear Solution (a = 0) 

V2 

1 

-l/x/2 

-1/V2 


Table 2.4: Summarized values from the similarity solutions. 


from which may be determined 


F(0) = 1.065... 

F„(0) = -0.881... 

F n (rhi P ) = -1.413... 

The result for F + (rj + ) is presented graphically in Figure 2.7. 

A summary of useful values for the problems considered above is provided in 
Table 2.4. 
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a 

b 

2a/ (1 + a) 

F+(0) 

.*£( 0 ) 

RSEM 

-1/3 

-1/3 

-1 

0 

OO 

OO 

-1/5 

-2/5 

-1/2 

1/4 

0 

2 

0 

-1/2 

0 

0.345 

-0.205 

1.702 

1/5 

-3/5 

1/3 

0.400 

-0.312 

1.581 

1 

-1 

1 

1/2 

-1/2 

V2 

a > 1 

-a/2 

2 

0.637 

-0.755 

1.253 


Table 2.5: Selected values from the generalized similarity equation, eqs 2.73 and 2.74. 


2.6.4 Generalized Similarity Solution 


It is instructive to digress here to point out several characteristics of the governing 
si mi larity ODE, eq 2.62. Restricting the discussion to cases a > — 1/3, selecting 
Ci = 1, C 2 = |5| 1//2 , and noting the condition b — — (1 + a)/2, the governing eq 2.62 
may be recast to give 


F + F+ +n 

77+TJ 


i+n+ 



2 a 

1 + CL 


■F + = 0 


(2.73) 


As described above, having introduced A, the boundary conditions can be specified 
at the tip, 


F + (l)~0 F+(l)==-l/2 (2.74) 

Solutions to eq 2.73 above for several values of a are plotted on Figure 2.8 and the 
significant numerical results are tabulated in Table 2.5. The cases of capillary rise, 
the spreading drop, and the constant flow condition are represented by the curves 
a = 0, -1/5, and 1/5, respectively, for this particular selection of C\ and €2- The 
solution F + = (1 — r? H_ )/2 is the asymptote for all solutions to eq 2.73 at the tip. 
Note that the leading order tip solution, F + = (1 — r] + )/ 2, found for the case a = 0 
(eq 2.66) is valid for any a . For a — l the linear solution holds, not only at the 
tip, but throughout r] + . Furthermore, the sole coefficient in eq 2.73-multiplying the 
last term of the left hand side-is only a slowly varying function of a: for a > 0 it 
is bounded between 0 and 2, and for all a the coefficient is bounded by -2 and 2. 
(Solutions for a < —1/3 produce curves which do not intersect the positive F + -axis.) 
Consequently, for a « 1 the solution F + (r} + ) is nearly linear. 
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Figure 2.8: Similarity function F + (r) + ) for a > —1/3. 

2.6.5 Inclusion of Gravity in ^-direction 

If gravity along the x-direction, g x , is incorporated into the analysis it can be shown 
that the governing partial differential equation is given by 

h T = 2hl(l H- 2 Bofjh?) + hh zz (l + Botfh?) (2.75) 

where Boh = fpgH 2 /c is the Bond number based on surface curvature (via /) and 
characteristic height H , and with ^-positive acting in the negative x-direction (see 
Appendix A. 4). By transforming this equation in like manner as for the similarity 
solutions above one obtains 

1— a— 26 

FF m {l+Bo H C 2 T 2 *F 2 )+2F${l+2Bo H C 2 T 2 *F 2 )--^-[aF + br,F n ] = 0 (2.76) 

which permits similarity solutions only when a — 0, 6 = —1/2. This means that only 
solutions of the type determined for the case of capillary rise are possible. For this 
choice the governing equation becomes 

FF m { 1 + Bo h F 2 ) + 2F 2 (1 + 2 Bo H F 2 ) + yF v = 0 (2.77) 

with h = F,ri = z{2t)~ 1/2 , and C = ??« p (2r) 1/2 , subject to boundary conditions 

F(0) = 1 F{r, tip ) = 0 
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and the integral mass balance 


rvtip _ 

2(1 + Bo h )F v (0) + / F 2 dr) = 0 

Since, in general, eq 2.75 assumes that surface tension effects dominate the flow 
(pgH 2 tan 2 a/ a <C 1), for the solution to be valid along the entire corner, it can be 
shown that the range of Boh is restricted by the condition f^Bon tan 2 a <C 1. In 
addition, the appropriateness of the boundary condition F(0) — 1 is unsubstantiated. 
Regardless, F is solved numerically for several values of Boh and the results are 
presented graphically in Figure 2.9. Solutions where g is negative (acting in the 
positive ^-direction) may also be determined and values for Boh > — 1 are shown. 
jF^(O) — * — oo as Bo h — > —1- Figure 2.9 illustrates the significant effect of negative 
g on 7^(0) and the relatively modest effect on r} tip . 7 The numerical values for both 
may be correlated with Boh to produce 

r) tip = 1.702 (1 + 0.286So//) 1/2 (2.78) 


F„(0) = -0.349 


(1+0.7Bo h ) 1/2 
1 + Boh 


(2.79) 


These relationships axe accurate to d=3% for the range Boh £ "1? but should 
be considered practical only in the limit f~ l Bo H tan 2 a <C 1. The case Bq h —>■ ► 0 
recovers the capillary rise result of Section 2.6.1. 


2.6.6 z-Curvature Consideration 


An interesting similarity solution is possible if the effects of curvature in the 21 - 
direction are considered in an asymptotic analysis. In the unique limits of e 2 <C 
e 2 / <C 1, it can be shown that 

Po = -<r (Jj + ( 2 fh Z zj (2.80) 

leading to the governing partial differential equation 

h T = 2 h 2 z H- hh zz - e 2 / (4 h 2 h z h zzz + h 3 h zzzz ^j (2.81) 

7 For the infinite column solution inclusive of g x , the fluid in the corner is unstable for Boh < — 1, 
see eq 2.58. 
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Figure 2.9: Similarity function F(rj) for the capillary rise problem for several values 
of Boh . 


Using the transforms of eq 2.59, similarity of the governing equation and boundary 
conditions may be achieved only by setting a = — b = 1 {h ~ rF, rj = z/r) and the 
resulting invariant system is 


. + 2F+ 2 +^-^ 


+ F + F*+ r)+r] + r]+ 


subject to the boundary conditions F + ( 1) = 0 and F++( 1) = —1/2. Two additional 
boundary conditions are needed and F+ +r] +(l) = F+ +v+rj+ (l) = 0 are likely good 
choices. In the limit e 2 / < 1, a closed form asymptotic solution is possible to the 
leading order system: choosing Ci — C 2 = 1 


F+F+ n+ + 2F+ 2 + 77 + F+ + - F + = 0 (2.83) 


subject to F + ( 1) = 0 and Fj+{ 1) = -1/2. A solution is 


F + = ^(1 - V) 


(2.84) 


from which F + ( 0) = 1/2 and ^+( 77 ) = -1/2. Using the integral mass balance 
equation Tf tiv may be determined, but only with knowledge of F( 0 ) which raises the 
question as to the nature of the condition of h( 0, r). If F(0) = 1, as was considered 
for the case of capillary rise, closed form expressions for h(z,r) and C{r) are readily 
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obtained. However, these solutions retain dependence on the length scales H and L 
and thus are not unique. A function T(e) is introduced such that .F(O) — T(e). The 
integral continuity equation thus gives 


Vu P = A = [r(e)] 1 ' 2 


-3 F f + 


7]+—0 


0 


1 , ' -! 
F +2 dv + 


1 1/4 


= [2r(e)] 1/2 (2.85) 


Solving for the dimensional values for h and L shows 


h = Gt 


fm 

2 


,. e [ r (*)] 1/2 

2 1 / 2 


( 2 . 86 ) 


where 


L = Gt- 


[r(e)P 

2*/2 


oFi sin 2 a 


(2.87) 

( 2 . 88 ) 


It is interesting that if F(e) = e -2 all H and L dependence vanishes, similarly as if 
e — 1. 

The leading order solution arises for the particular selection a — — 5 = 1. This 
choice was necessitated by the presence of the 2 -curvature terms, which being 0(e 2 /), 
are eventually ignored. The solution is not unique as it satisfies the homogeneous 
system of eq 2.37. This linear result is shown on Figure 2,8 (with a = 1). 


2,6.7 Summary of Closed Form Results (dimensional) 


The analytical and numerical solutions for the cases of Sections 2.6.1-2.6.5 yield 
useful constants which allow the closed form determination of particular quantities 
of interest. These include primarily tip location and velocity (by differentiation), but 
flow rate and meniscus height and velocity at 2 = 0 may also be determined if desired. 
The following are some of these results derived from the similarity solutions: C is the 
tip location, h is the meniscus height and Q is flow rate at z = 0. All quantities are 
dimensional in this section and the coefficient G is given by eq 2.88. Fa is given by 
eq 2.35 and the steady solution of eq 2.52 is provided for completeness. (Rapid and 
accurate estimations of the quantities listed may be made by noting Fa ~ tan a and 
Fi « 1/7.) 
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Steady Solution, h = h(z). For the special case of h( 0) = Hi, h(£) = (C = 
const), eq 2.52 produces the dimensional results 


h = H i 



(2.89) 


HI)’ 


(2.90) 


the latter which gives a ‘maximum’ capillary flow rate along a corner for these con- 
ditions. 


Capillary Rise, h( 0, t) = H. The closed form results for this problem contain H which 
is a quantity assumed constant for the flow. The determination of H for a number 
of simple test cell types is discussed in Section 2.6.8. 

£ = 1.702 G 1/2 H i/2 t l/2 (2.91) 

Q = 0M9F A G 1/2 H 5/2 r 1/2 (2.92) 

If the effect of g along the x-axis is considered, the above are modified by the corre- 

lated coefficients for r\ tip and F v (0) given in eq’s 2.78 and 2.79, which when substituted 
produce 

C = 1.702(1 + 0.286 Bo H ) 1/2 G 1/2 H 1,2 t l/2 (2.93) 

Q = 0.349 ^ + F A G 1/2 H 5/2 t~ 1/2 (2.94) 

1 + Boh 

limited to Boh ~ — 1. 


Spreading Drop, Volume H 3 . Since an exact solution exists for this case, h itself can 
be expressed in closed form. 

h = 0.70 6Ff /5 G~ 1/5 H 6/ h- 1/5 (l - v 2 />?) (2-95) 

where A = (15/2) 1 / 5 and 0 < r\ < A with 

r) = 0.796 F\ ,5 G- 2/5 H~ 3/5 r 2/5 z 
C = 1.87 9i^ 1/5 G 2/5 i7 3/5 t 2/s 
h{0,t) = 0.706FJ 2/5 G- 1/5 # 6/5 r 1/5 


(2.96) 

(2.97) 


Constant Flow Q imposed at z = 0. 


C = 1.634 (jFj 1/5 G 2, H 3,z 


h(0,t ) = 1.279 


'A) 

Fa) 


2/5 


G -l/5 t l/5 


(2.98) 

(2.99) 
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2.6.8 Determination of if, Capillary Rise 


As will be shown from experimental results, there exists a location zjj such that the 
interface height h(zn , t) appears to rapidly rise to and remain fixed at a certain height, 
H. It is argued by Dong and Chatzis [17] that a constant pressure condition is estab- 
lished in the fluid at a certain ^-location for large time. Prom this constant pressure 
condition the constant height H may be evaluated. Using the elegant approach of de 
Lazzer and Langbein [29] , H is determined below for the n-sided regular-polygonal 
and rectangular cross-sectioned cylinders with sharp corners of face widths D and 
D i) respectively: 


Hn—poly — 


D sin(a 4- 8) 

V F An 


y /2 

\ sm 2 (a4<5 ) ) 


( 2 . 100 ) 


where 8 — n/n — 0, a = n (n — 2)/2n, and 

rr jD 2S m(a + S) [, D 1 4 F An \ 1/2 ' 

4/ F An V D 2 fsin 2 ( a + 6)) 


( 2 . 101 ) 


Note for the rectangle that D\ < D 2 , j = 1 4 D\jD 2 , and H n - po i y = H rect for the 
square cross-section, n — 4, Di/D 2 = L Fa u is given by 


_ sm 2 8 . c c „ 

FA n = - — — 4 sm 6 cos 8 — 8 
tana 


( 2 . 102 ) 


For the limits of 0 < Di/D 2 < 1 it can be shown that 1/2 < H rect f cos 0/Di < 1. 
Since ultimately C a if 1 / 2 , the effect of aspect ratio on C is limited to a ~ 25% 
increase over that of the square cell. 


The scaling introduced in Section 2.1 restricts Fi to a narrow range of values 
1/8 ;$ F{ < 1/6. Thus, G ~ sin 2 a/f for all values of a and 8. The dependence of 
the tip location C on system versus interface geometry can now be discerned. Recall 
that for this problem of capillary rise the tip location is £ ~ (GH) 1 / 2 for the n-sided 
polygon. For 

F A n cot 7r/n 


sin 2 (a 4 8) 


<1 


(2.103) 


it can be shown that H ~ 1/f sin(a 4 8) for a fixed. For a fixed it can also be shown 
that G ~ 1//. Thus L ~ l/f where the contribution from H 1 ^ 2 ~ / -1 / 2 is strictly 
geometric and is tied to the cross-flow area whereas the contribution G 1//2 ~ f~ 1 ' 2 
stems from the driving pressure gradient. In this way the effect of increasing contact 
angle, or equivalently the limit <5 — > 0, is seen to impact the tip location C by reducing 
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the driving pressure gradient as well as the cross-sectional area for flow. This is also 
true for the flow rate Q since Fa is only a weak function of 6 for a fixed. 

The reason the cross-sectional area is affected to this degree by changes in 8 is due 
to the system pressure argument leading to the constant height boundary condition 
for this problem. For flow scenarios other than the capillary rise problem the effects of 
increasing contact angle are predominantly due to reductions in the driving pressure 
gradient and not to changes in cross-sectional area. This should be evident by the 
fact that Fa is a banded function of <5, Fa ~ tana for all <5, and that A' = H' 2 Fa ■ 
Therefore A 1 exhibits weak dependence on S. As a case in point, for the problem 
of the spreading drop (Section 2.6.2), H is the cube root of the drop volume and is 
unrelated to any system parameters. The effects of increasing contact angle for this 
problem come primarily through / as found in G (see eq 2.95) since Fa is a banded 
function in terms of 8. 


2.7 Quantified Constraints 

The solutions obtained above rely on a rather large number of assumptions, the 
validity of which will now be investigated. The solutions themselves will be used 
to estimate the order of magnitude of the terms neglected in the analysis. A list of 
assumptions for the capillary rise problem is followed by a similar list for the problems 
of the spreading drop and a constant flow. Each assumption is expressed using the 
dimensionless numbers of eq 2.104 when possible. Otherwise dimensional quantities 
are used. 

While the asymptotic analysis formally requires only e 2 — * 0 and e 2 / 1, for 

practical application it is prudent to consider the order of magnitude of the other 
parameters which might appear in the problem. Thus, the constraints on 7£, gravity 
level, and the contact line boundary condition are also considered below. 

For these calculations the following terms will be used to simplify the expressions: 

Boh = Ca t = ^r T=—G- (2.104) 

a crt Fi sm a 

where gi refers to the component of g in the i-direction. The lead constants ap- 
pearing in the results are known exactly and are rounded and provided in decimal 
form for completeness. It is important to note that for moderate a in the range 
7r/12 < a < 57r/12, and for 8 >, 7r/18 (implying lower contact angles), T (oc T) is 
approximately an 0(10)-0(100) quantity (see also Figure 2.6). As 8 — ► 0, T — ► oo. 
On the other hand, / is ~ 0(1) for moderate a and 8 >, 7r/18, but tends to ~ 0(0" 2 ) 
and 0(a/ sin 8) as a tends to n/2 and 0, respectively. 
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The standard assumptions of a passive overlaying gas, constant properties, and 
ideally sharp interior corners of infinite extent are not addressed analytically. 


2.7.1 Assumptions for Capillary Rise 


• e 2 < 1 [e — H/£(t)]. This is the slender column assumption and requires 

OMSCatF < 1 (2.105) 

This requirement is not appropriate near the bulk meniscus where e ~ 0(1), 
locally. 

• e 2 / < 1. This is the requirement of slight surface curvature along the corner 
axis as compared to curvature perpendicular to it. It is given by 


OMbfCatF < 1 (2.106) 

This condition is difficult to satisfy as 8 — ► 0 since / — > oo in this limit. T 
behaves similarly. 


• 1Z 1, or rather 7 Z ~ 0(e 2 ) (see eq 2.7). In effect, this constraint is necessary 
in order to neglect inertial effects on the left hand side of eq 2.4. For prac- 
tical low-gravity systems of typical dimensions and fluid properties it usually 
cannot be satisfied. It is often difficult to satisfy for normal-gravity systems as 
well. A more precise accounting of the constraint of low inertia is offered by 
\lZw 0t /P z \ <C 1, a quantitative comparison of the unsteady term to the driving 
pressure gradient. The convective term is found to be an order of magnitude 
less than the unsteady term (\w o w Oz /w 0t \ ~ 0,1) and so can be ignored. This 
constraint can be written as 


k 2 sin 2 a 
2vt 


<1 


(2.107) 


A worst case is when h = H. The presence of inertia in the problem will result 
in the retardation of the flow, particularly during start-up. This constraint is 
quite limiting for applications in low-gravity environments where H ~ O(lOcra) 
might be anticipated. However, near the tip where h — > 0, inertia is rightfully 
ignored save at very small times. 


g x ~ 0. Negligible gravity in the x-direction requires 


(2.108) 
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• g z ~ 0. Gravitational effects parallel to the corner axis (^-direction) are negli- 
gible when 

Bo hc = = 1.702 Bo h {Ca t P)~ 112 < 1 (2.109) 

C 7 

The impact of g z increases with time as ~ t 1 ^ 2 . 

• h{ 0, t) = 1, constant height assumed at z = 0 for all t. This initial condition 
was selected based on experimental observations. Its validity needs be verified 
experimentally and will be addressed in Chapter 3. 

• 0 = const. The solution for w 0 may be decomposed into portions perpendicular 
and parallel to the contact line. The component perpendicular to the contact 
line may be used in turn to determine the local Capillary number {Cad), the 
magnitude of which may be used to estimate departures from the constant, 
static contact angle assumption. Capillary numbers of sufficient size signal 
the dynamic wetting phenomena leading to a dynamic apparent contact angle 
{9 d ) effect which can decrease (increase) the effective curvature (contact angle) 
of the free surface and so reduce the pressure gradient driving the flow. For 
advancing and receding contact lines, the correlation of Hoffman [30] may be 
used 

6 d = 4.54(5 + Ca d )°' 353 (radians) (2.110) 


where 


5 = 



( 2 . 111 ) 


provided B + Ca d > 0. The constant B incorporates the effects of partial 
wetting. Firstly, Ca d is determined using 0 = 6 stat . This value of Ca d is then 
substituted into eq 2.110 and 6 d is estimated. Solutions producing significant 
changes between 6 $ta t and 0 d thus determined should be reconsidered in light 
of dynamic contact angle effects. 

The dimensional velocity of the contact line V d is determined from 

v d = (h t -(w 0 )h z ){i+hiy 1/2 


Using this form and solving for Ca d = (jV d /a gives 

Ca d = h 2 n Ca t (l - (l + (2.112) 

Since h for this problem was determined numerically, h v must be evaluated. 
However, h v is bounded such that -0.851 < h v < -0.349. If the constraint 
e 2 <C 1 is satisfied using eq 2.105, and provided T £ 0(1), 8 eq 2.112 solved at 

8 This will likely satisfy the constraint fe 2 < 1 as well. 
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the tip (the worst case) gives 


Ca d = 1.44 8Ca t (2.113) 

which can be used to determine 8 d from eq 2.110. 

A quick and worst case estimation of dynamic contact angle effects may be 
made on the tip velocity prediction using 


VtiplOstat) - (FA (f\ 

v tip (e d ) uA stot U/ 



(2.114) 


which gives the upper limit in terms of percent error due to 0 d . The error 
decreases with increasing t through the effect of Ca t oc 1/t on Ca d and thus 
8 d . It is also worthy of note that for fixed a, / is more sensitive to changes in 
8 d than is F x . For instance, for a = 7r/4 a change of 8 d from 0 to 7 r/6 yields 
an increase in Fi of 10% while f increases 85%. As 8 — ► 0, / — > oo, while F* 
cannot fluctuate more than 25% (1/8 £ Fi < 1/6). Therefore the principal 
influence of 0 d is not through Fi, but through / originating in the pressure 
scale, P' ~ cr/h'f. 


2.7.2 Assumptions for Spreading Drop 


For these calculations Ca t , T, and Boh are the same as defined above and Ii is 
determined from the drop volume (i? 3 ). In addition, A == r} tip and is given by eq 2.69 
as determined by the exact solution to the similarity equation. Thus, 0 < rj < A. 
Since Fa ~ tana, in the limits a — > tt/2, or a — ♦ 0, then Fa ~ (tt/ 2 — a) -1 , or 
Fa ~ a, respectively (see Figure 2.3). 


e 2 « 1 [e = h(z, t)/C(t)\ is the slender column assumption. 

, 2\2 


0.141F/ /5 (C'a t .F) 6/5 ^1 - ^ < 1 


(2.115) 


The worst case is when rj = 0 

• e 2 f 1 requires slight curvature along corner axis compared to perpendicular 
to it, 

/ 2 \ 2 

O.UlfF^iCatF) 6 ' 5 ( 1 - ) <1 (2.116) 
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1Z ~ 0(e 2 ), or \TZw 0 JP z \ <C 1 is the condition necessary to neglect inertial 
terms. 



\vt) f\ \ \y 


(2.117) 


This term is zero at the tip ( r} tip = A) and is maximum at g = 0 where inertial 
effects are greatest. 

g x ~ 0 is the assumption that the effect of gravity perpendicular to the corner 
axis (x-direction) is small. This constraint can be written as 

Bo h = = 0.5Bo h Fa 4,s {Ca t Ff 5 (1 - <1 (2.118) 

The effect of g x diminishes as the tip is approached. 

g z ~ 0 is the assumption that the effect of gravity parallel to the corner axis 
(^-direction) must be small, 

Bo hc = = 3.0 Bo H F^ 3/b (Ca t F)~ 1/5 ^1 - ^ < 1 (2.119) 

The effects of g z decrease as the tip is approached, but increase weakly with 
time as ~ t 1 ^. 


0 = const. As for the ease of capillary rise, departures from 6 stat may be esti- 
mated by determining Ca d> The steps follow similarly those outlined for the 
capillary rise problem, only here 


Ca d = 0.42FJ 2/6 (Ca^) 1/5 ^ - lj 




- 1/2 


( 2 . 120 ) 


It is interesting to note the location where Cad = 0, namely, 77 = A/3. The 
contact line upstream of this location is receding while that downstream is 
advancing. When 0.16 (Ca t P) 2 <C 1 the above reduces to 

Ca d = 0A2Ff /5 Ca t ( Ca t F ) l/ 5 - 


( 2 . 121 ) 


Hoffman’s correlation may be used to estimate 64. 



2.7.3 Assumptions for Constant Flow 


For this section the introduction of the length scale 


(ixTtQ 2 

V ° F a 


1/5 


( 2 . 122 ) 


simplifies the expressions and is used as l H' was above in Ca f , Bo^-, etc, This 
value for H is time dependent and increases with t 1 ^, Note that H increases with 
increasing viscosity and flow rate Q. 


• e 2 <C 1 [e = h(z,t)/£(t)] is the slender column assumption. The worst case at 


V = 0, 


0.61 



< 1 


(2.123) 


Although H increases with time, the length of the column increases faster and 
e 2 oc r 4 / 5 . 


• e 2 /<l describes slight curvature along corner axis compared to that perpen- 
dicular to it, 


0.61/ 



«1 


(2.124) 


® 7 Z ~ 0(e 2 ) or \TZw 0 JP z \ <C 1 is the condition necessary to neglect inertial 
terms. It is zero at the tip and maximum at rj = 0. At rj — 0, 


2.62 


sin 2 a H 2 
/ vt 


<1 


(2.125) 


Thus the constraint decays slowly with time as t 3 ^ 5 due to the increasing 
height H . 


• g x ~ 0 describes the assumption that the effect of gravity perpendicular to the 
corner axis (x-direction) must be small. This is always true at the tip but for 

77 = 0 , 

Bo h = = iMBoh < 1 (2.126) 

a 

Note that Boh increases with t 2 ^ and Q 4/5 . 
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g z ~ 0 gives the assumption that the effect of gravity parallel to the corner axis 
(^-direction) must be small and satisfy 


B0hc = = 2.09 Bo h * 2/5 « 1 (2.127) 

which shows Bohc increases with £ 4 / 5 and Q 3 / 5 . 

• 0 = const. The presence of 0d effects may be estimated similarly as for the 
previous cases. For the entire domain, 

Ca d = 0 A8FiCatF^ (l - ^ [l + 0A8Ca t FFj] ~ 1/2 (2.128) 

At the tip 

Ca d = 1.92 FiCat [1 + 0.958C'a i .7 r r 1/2 (2.129) 

and for small Ca t J - with Fi cz 1/7, the above reduces to 

Ca d = 0.27Ca t (2.130) 

Ca d oc £ _4 / 5 (/i/cr) 6 / 5 Q 2 / 5 and may be used with Hoffman’s correlation to es- 
timate 0 d . The capillary number for this problem is higher in magnitude and 
decreases at a slower rate than for the problems of capillary rise and the spread- 
ing drop. 

2.7.4 Abbreviated Summary of Assumptions 

The assumptions analyzed above are summarized in Table 2.6. In this table, values 
for Tj are chosen at the origin or at the tip ( rj = A) in order to present the maximum 
magnitude of the constraint under consideration. The lead constants are dropped for 
clarity. It should be noted that the condition to neglect inertia is oc 1 j\x while that 
to satisfy e 2 1 is oc \x. 
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Constraint 

Capillary Rise 

Spreading Drop 

Constant Flow 

e 2 < 1 

Ca t T 1 

FX Vi (Ca t J r ) e,s < 1 

w « i 

e 2 /< 1 

fCa t T<C 1 

fFf /5 (Ca tFf' 5 < 1 



if 2 sin 2 a ^ i 
i/t ^ 

1 ' (dH*\ 2 ( H 2 \ 3 Ff sin 3 a 

[Ov UJ -f=^r\ 

1/5 

<1 

if 2 sin 2 a ^ i 
/vt ^ 1 

~ 0 

Bo h -C 1 

Bo h FX 4/ 5 (Ca t Ff /5 < 1 

Bo# < 1 

9z ~ o 

see eq 2.119 

Bo h F^ 3/ - (Ca t F)- 1/b < 1 

see eq 2,127 

h(0,t) = 1 

see Chap. 3 

- 


- 

^(C'acz) for 
e 2 < 1 
Odj <C! 1 

Cad ~ Ca t 

Ca d ~FX 2,5 Ca t {CatF) 1/!i 

Cad rv ' Cuj 


Table 2.6: Summary of constraints. 
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Chapter 3 

Capillary Rise Experiments 


In this chapter a brief introduction to experimentation in a low-gravity environment 
is followed by a description of the capillary rise experiments and test procedures. 
A summary of a small portion of the test results is presented in graphical form to 
illustrate several global features of the flow. Digitized surface traces of the devel- 
oping flow are compared with the analytical results of the preceding chapter. The 
initial transient and the ‘start-up’ flow for a sample test is discussed in light of ex- 
perimental observations. A simple method for estimating tip locations at small times 
is suggested. For the large capillary systems which can arise in low-g environments 
one might expect significant inertial effects. However, it is observed that viscous 
forces dominate the flow in the tip region to a high degree, and that the tip flow is 
somewhat independent of a container length scale at small times. 

Extensive experiments for the capillary rise problem were conducted which quan- 
tify the broader effects of corner angle, container size, interior surface conditions, 
and a wide range of fluid properties including system contact angle. These will be 
described in a subsequent publication. The objective for presenting only a sample 
of these results here is to illustrate the utility of the analytical solutions. A host of 
normal and low-gravity tests were also performed for the case of a spreading drop 
of fixed volume in an open corner (see Section 2.6.2). These tests were performed 
to note the impact of corner angle, contact angle, drop volume, gravity, and fluid 
properties on the capillary driven flow. These results will also be described in a 
subsequent publication. 
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3.1 Introduction to Low-g Experimentation 


Experimentation in a low-gravity (low-g) environment is unique in several respects. 
Container sizes can be orders of magnitude larger than their normal-gravity coun- 
terparts. The large cell size, in addition to extending the range of certain fluid 
parameters, also allows greater precision in container fabrication, diminishing the 
role of surface contamination through corner irregularity and roundedness, surface 
roughness, and other heterogeneities. In addition, larger scale capillary flows are 
not difficult to visualize and can be tracked over long distances. The near absence 
of gravity reduces the need to consider its complicating effects in analysis. Several 
important applications of capillary driven flows exist in the low-g environment, e.g. 
fluids management aboard orbiting spacecraft. It is appropriate to study such flows 
in a low-g environment. 

A 2,2 second drop tower test facility at NASA’s Lewis Research Center is used 
to access the low-g environment for the capillary rise tests. Though the low-g time 
afforded by the drop tower is limited, the cells tested are sized to make adequate 
use of the time available. The 2.2s drop tower is approximately 27m in height. To 
minimize the ill-effects of aerodynamic drag on the experiment during free fall a drag 
shield (a box-in-a-box) technique is employed. Gravity levels of 10 ~ 5 g 0 are common. 
A detailed description of this facility may be found in [31]. The description here will 
be restricted to the test/drop rig apparatus. 


3.2 Description of the Experiment 

A transparent test cell is first mounted to a test stand and installed in a drop rig 
apparatus. Prior to release of the apparatus into free fall, a prescribed amount of 
fluid is injected into the cell, partly filling it. Upon release, hydrostatic forces are es- 
sentially eliminated and capillary driven flow results. The flow process is backlighted 
by a diffuse light source and filmed at a long working distance with a high speed 
cine camera at 128 frame/sec. A sample of frames taken from the film records is 
provided as Figure 3.1 for two tests of silicone oil fluid (PDMS) in a 22.6mm equilat- 
eral triangular cell. Figure 3.1a is for 2.0cs fluid at 0.2s intervals while Figure 3.1b 
is for lO.Ocs fluid at 0.5s intervals (the top cross-sectional view of these cells is that 
of Figure 3.2b). Quantitative data is digitized directly from the film records. Up to 
four test cells could be tested simultaneously. 

The principle variables for the drop tests are the container face width D and fluid 
viscosity fi. Corner half-angle a and liquid-solid contact angle 9 are also parameters 
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Figure 3.1: Selected frames from cinefilm records of capillary rise in equilateral tri- 
angles, D — 22.6mm. (a) 2.0cs fluid at 0.2s intervals (Tl), (b) lO.Ocs fluid at 0.5s 
intervals. The top view of the cells is that of Figure 3.2a. The flow along the right 
corner is analyzed. 



Fluid 

v ±2% 

(cs) 

p ±5 
(kg/m 3 ) 

<7 ±5% 
(N/m) 

No 

PDMS 

0.65 

760 

0.0159 

1.374 

PDMS 

1.0 

816 

0.0174 

1.383 

PDMS 

2.0 

872 

0.0187 

1.390 

PDMS 

5.0 

913 

0.0197 

1.396 

PDMS 

10.0 

935 

0.0201 

1.399 

PDMS 

20.0 

949 

0.0206 

1.401 


Table 3.1: Test fluids and properties. Polydimethylsiloxane (PDMS, silicone oil) 
provided by Dow Corning. 0 stat = 0 on all surfaces. 


Cell Cross Section 
(section type) 

Material 

Nr> 

(5893A) 

L 

(mm) 

D or Di x D 2 
(±0.05 mm) 

a 

(deg.) 

equil. triangle (b) 

acrylic 

1.491 

150 

12.0 

30 

equil. triangle (b) 

acrylic 

1.491 

150 

22.6 

30 

near square (c) 

acrylic 

1.491 

150 

6.1 x 6.7 

45 

near square (c) 

acrylic 

1.491 

150 

12.9 x 13.2 

45 

rectangle (c) 

acrylic 

1.491 

150 

6.0 x 13.0 

45 

rectangle (c) 

acrylic 

1.491 

#0 

8.1 x 25.6 

45 

equil. triangle (b) 

quartz 

1.459 

60 

5.0 

30 

equil. triangle (b) 

quartz 

1.459 

60 

10.0 

30 

square (a) 

quartz 

1.459 

60 

5.0 

45 

square (a) 

quartz 

1.459 

60 

10.0 

45 


Table 3.2: Test cell data for capillary rise tests. 


of interest, but in this report only systems with 9 $ta t = 0 will be considered. The 
properties of the liquids and test cells appropriate for the tests to be presented are 
listed in Tables 3.1 and 3.2, respectively. 

As the physical condition of the corner is critical to such flows the manufacture of 
the test cells is also critical. Figure 3.2 displays the joining technique of the faces for 
the cells. The acrylic cells are glued and the quartz cells are fused at the joints. In 
order to assess fabrication variability, four reproductions of each quartz section type 
and two of each acrylic cell were constructed. Tolerances of the critical dimensions 
are included in Table 3.2. Variations in corner angle are minuscule. 

The dim ensions for the rectangular cells are also included in Table 3.2 where 
L here is the total length of the cell. The majority of the tests were performed 
using cells of square or equilateral triangular cross sections. For both of these cell 
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Figure 3.2: Cross sections of capillary rise test cells. 


configurations, the cells are fabricated such that the meniscus centerline along the 
corner h(z : t) can be viewed in profile with a minimum of optical distortions. Figures 
3.2a and 3.2b show the cross-sections for these cell types. A ray trace analysis, 
described in Appendix B, which provides corrections for distortion when refractive 
indices are not matched, is employed for certain measurements. These corrections 
are never greater than 10%. 

Cell cleaning procedures were developed for each cell material to provide an ac- 
ceptable degree of control over the surface condition. The cell preparation includes 
washing with a strong soap solution in hot water, rinsing thoroughly with distilled 
water, then rinsing with ethanol and/or methanol followed by a final distilled water 
rinse. The cells are allowed either to air dry or are dried in a vacuum oven. An 
additional acetone wash step is added for the glass cells before the ethanol rinses. 


3.3 Preliminary Measurements 

The principal measured quantities from the film records are the location of the lead- 
ing tip of the fluid in the corner £(t), the global meniscus height Z(t), and the surface 
profile, or meniscus centerline height h(z, t) (see Figure 3.1 for notation, all length 
measurements are in mm). A sample of these results is provided and, where appli- 
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cable, the dependent variable is plotted against t 1 / 2 as suggested by theory. Figures 
3. 3-3. 7 are intended to convey the general character of noteworthy aspects of the flow 
using the point location measurements of £ and Z. A comprehensive analysis of two 
data sets is deferred to Section 3.4. 


3.3.1 Tip Location, C(t) 

Figures 3.3 and 3.4 display data comparing the impact of viscosity, cell size, and 
corner angle on the tip location, £(t). Figure 3.5 illustrates the effect of container 
aspect ratio for rectangular cross-sectioned cells. All tests here are performed with 
8 stat — 0 in acrylic cells. The linearity of the data seems to imply the t 1 / 2 dependence 
of £ for this flow which was demonstrated by Dong and Chatzis [17] to be correct at 
long times. As anticipated, the tip velocity increases with decreasing viscosity and 
increasing cell size as observed in [17]. The tip velocity is also seen to increase for 
a = 30° compared to a = 45° for cells of approximately the same inscribed circle 
diameter. 1 The solid lines are fits £ =» c + rat 1 / 2 to the data (symbols) and exclude 
the initial data points for times t ;$ 0.1s. 

Concerning Figure 3.5, it is interesting to observe that the slope of the fits in- 
crease only modestly for a 4-fold increase in aspect ratio. This is due to the relative 
insensitivity of H to Di/D 2 as could be demonstrated analytically by eq 2.101. 

The measurements for these figures are made using a film motion analyzer in 
manual mode without image processing. Multiple drop tests are made to note the 
repeatability of the measurements in regards to cleaning technique, fill technique, 
etc . The majority of tests show repeatability 2 to within ±6% when comparing fitting 
parameter coefficients for the different runs and ±3% is common. With this level 
of precision it is easy to quantify differences in flow behavior from vessel to vessel, 
and even from corner to corner of the same vessel. Concerning the latter, as regards 
the tip location, asymmetric flows are detected between different corners of the same 
vessel. These differences are repeatable (again to within ±6%) and are linked to 
discrepancies in the fabrication quality of the corners. The largest differences are 
observed in the tip location measurements for the rectangular cross section vessels, 
yet are always less than 10%. 

It is important to note that when the Concus-Finn condition is satisfied, regardless 
of the gravity level, a thin column of fluid already exists in the corner prior to the 
drop test [2]. This preexisting film makes it difficult to determine the true tip location 

1 Recall that the case of a = 30° provides the most rapid response of the fluid (see Section 2.5). 

2 This is true when the number of duplicate tests was statistically significant. 
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Figure 3.3: £ (mm) as a function of t 1/2 (s 1/2 ) in acrylic triangle, D = 12.0mm, for 
fluids of different viscosity. □ 1.0, A 2.0, Q 5.0, and V lO.Ocs. 
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Figure 3.4: £ (mm) as a function of t 1 / 2 (s 1 / 2 ) for a fluid of v — 5cs in an acrylic 
cell showing the effect of cell size and corner angle. □ D = 22.6mm triangle, Q 12.0 
triangle, A 12.9 near square, and V 6.1 near square. 
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Figure 3.5: £ (mm) as a function of t 1/2 (s l/2 ) in acrylic rectangular cells for a fluid of 
v — 5.0cs showing the effect of aspect ratio. Q Di*D 2 = 6.1 x 6.7mm, V 6.0 x 13.0, 
and □ 6.1 x 25.6. 


of the meniscus during the early stage of the drop tests. An apparent C{t) was thus 
measured at the onset of the flow by extrapolating the apparently linear meniscus 
profile near the tip to the vertex of the corner. The infinite rise of the fluid in the 
corner is based on a continuum approach and of course an infinite corner. Provided 
the corner is long enough, the effects of corner imperfection, finite molecular size, 
Van der Waals forces, and even evaporation will prevent such an unbounded height. 
As the fluid column rises during the drop, upper portions of the corners are reached 
where the film is not observed and data are taken in these regions without any need 
for extrapolation. A smooth transition between data measurement techniques is 
achieved as can be seen from the data. It is arguable that flow in the preexisting film 
is so dominated by viscosity, and the slope of the surface so slight, that its presence 
for a portion of the bulk flow (on average 0.1s of the 2.2s drop time) is of little 
significance. 

Another complication in the tip measurements is that during the latter stages of 
the corner flow the tip location for some of the tests (usually low viscosity fluids in the 
smaller test cells) becomes increasingly difficult to locate as the slope of the meniscus 
approaches zero. Though the rise can be viewed without mistake at high playback 
speeds, stop-action digitization of the tip location leads to some loss of accuracy for 
these cases. How much loss is uncertain. For these tests, repeat measurements of a 
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Figure 3.6: Z (mm) as a function of t 1 ^ 2 (s 1 / 2 ) in an acrylic triangle with D — 12.0mm 
for fluids of various viscosity. O ^ = 6-65cs, V 2.0, □ 5.0, and A 10.0. 

single test event show measurement repeatability (‘experimenter’ error) on the order 
of the experiment repeatability itself, namely ±3-6% when comparing fit coefficients. 

3.3.2 Global Meniscus Height, Z(t ) 

Figures 3.6 and 3.7 show the decrease of Z with time, which is related to the drain 
rate of the container. As was the case in Figures 3.3 and 3.4, for the cells tested here, 
the flow rates and flow velocities increase with decreasing viscosity, increasing cell 
size [17], and decreasing corner angle for this range of a >, 30°. The dimensions of 
the cells used in Figures 3.4 and 3.7 were selected so that the radius of the inscribed 
circle is approximately the same for the large triangle and the large square cell. The 
small triangle and small square cells also share approximately the same inscribed 
circle radius. 


3.4 Detailed Measurements 

From the results of the point measurements presented in Figures 3.3-3. 7 one might 
conclude that the similarity solution of Section 2.6.1, producing the t 1 ^ 2 dependence 
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Figure 3.8: £ (mm) as a function of t (s) in acrylic triangular cells with c+m(t-t f ) : 1/2 
fits, v v — 2.0cs, D = 22.6mm (Tl: c = —6.17mm, m — 54.88mm/s 1 / 2 , tf = 0.041s), 
O v = 10.0, D = 12.0 (T2: c = -0.43mm, m = t f = 0.063s). 


for the flow, is in general correct for times t £ 0.1s. Such a conclusion is premature 
as will demonstrated below using surface profile measurements of h(z, t). 

The measurements from two capillary rise tests of PDMS in equilateral triangular 
cells will be presented and contrasted in this and the following sections. In the test 
denoted Tl the viscosity is 2.0cs and the face width D = 22.6mm. For the test 
denoted T2 the viscosity is lOcs and D = 12.0mm. C(t) and Z(i) are presented in 
Figures 3.8 and 3.9 for both tests, the data of the former being fitted by a line of the 
form C — c + m(t — t/) 1 ^ 2 . 


3.4.1 Surface Profile, h(z,t) for data set Tl 

In Figure 3.10 are overlayed ten digitized surface traces for the corner flow of test 
Tl. These data, and the data to be presented subsequently, are acquired using an 
automated image analysis system which digitizes the cine film records for selected 
frames and tracks the development of the meniscus in time. The grabbed images are 
low-pass filtered and a median threshold is selected to define the interface which is 
then skeletonized. Figure 3.10 suggests that there is a location z = z H such that 
the height at the location h(z H , t) = H is fixed for sufficiently long times. In the 
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Figure 3,9: Z (mm) as a function of t (s) for T1 (v, t 0 — 0.09s), T2 (□, t Q ~ 0.03s). 


figure, this location appears to occur at approximately z H ~ 11mm. A magnification 
of the region near this point depicts the waist created by the curves passing through 
a narrow neighborhood near Zh, Figure 3.10b. Measurements of the local surface 
height, h(zi,t), at several ^-locations, Zi, in this region are presented in Figure 3.11. 
The distance between each Zi is approximately 1.5mm. As can be seen from the 
figure, depending on z i) h{z^t) will in general under- or overshoot h(zn,t) = H. 
The data of Figure 3.11 for t > 0.7s are replotted and fit with linear coefficients 
on Figure 3.12. This figure reveals that there indeed exists a z = zh such that 
h(z H ,t) — H exp = const (H exp = x wa ii — %h)- Interpolation between the smallest 
positive and negative slopes provides the empirically determined z — zh and H exp . 
The interpolated zh is between the hollow and solid triangles on the figure. In 
subsequent figures, the origin is moved to z = zh • This location then defines the 
apparent origin where the constant height boundary condition of the asymptotic 
analysis can be applied. Downstream of this location h is always increasing towards 
H while upstream h is decreasing everywhere to this same level. Since only limited 
time was available for these tests the long time behavior of the boundary condition 
at Zh cannot truly be confirmed. However, the results of Dong and Chatzis whose 
experiments lasted hours offer strong evidence that the constant height (‘constant 
pressure 5 ) condition is acceptable. H = 3.67mm determined analytically as discussed 
in Section 2.6.8 compares well with H exp = 3.73 ± 0.15mm determined experimentally 
for this test, as noted on Figure 3.11. Uncertainty in the experimental measurement 
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30 35 X 40 38 39 X 40 


Figure 3.10: Traces of the surface height h(z,t) measured along x for Tl. As time 
increases the profiles recede at the centerline x ~ 30mm and advance along the wall 
at x — 43.0mm. The surface profiles are shown for times t = 0, 0.078, 0.156, 0.234, 
0.312, 0.390, 0.547, 0.703, 0.859, and 1.015s. (a) Global profile to cell centerline, (b) 
enlarged view of dashed box area on (a). The coordinates x and z are in mm. 
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0.00 0.25 0.50 t 0.75 1.00 

Figure 3.11: h(zi,t) (mm) as a function of i (s) for Tl. Each data set is taken for 
constant Zi with an increment of approximately 1.5mm between each Z{. 

of H is due to the optical scale factor of the measurements and the corrections for 
image distortion caused by refraction. 

H is not fully established until approximately in — 0.4s. The initial inertial tran- 
sient t 0 ~ (pRl/a) t/ * for the establishment of the low-g interface may be observed 
by the rapid decline of the meniscus Z{t) early during the drop test. These times 
are noted on Figure 3.9 along with the estimates t Q — 0.09 and t 0 — 0.03s for Tl 
and T2, respectively. Note that t Q ~ 0.7b(pR 3 Jo) 1/2 and that the general order of 
t 0 agrees with the fitted time origin tf of Figure 3.8. Siegert et al. [32] demonstrate 
that the initial time response of an interface after a step reduction in gravity level 
may be correlated by t 0 = C(pR 3 /a)^ 2 , where R is a characteristic length and C 
is a constant, both of which depend on the container geometry. For a right circular 
cylinder of radius R — Ri C they find that C ~ 0.413. 

The experimental values of zh and tn are used to rescale the height data; Figure 
3.13 shows h as a function of zt~ x ! 2 and z(t — t#) -1 / 2 where z — Zi~ z H . The failure 
of the curves to collapse when plotted against the ( dimensional) similarity parameter 
zt~ l / 2 suggested by theory attests to the significant effects of the start-up transients 
necessary to establish the constant height condition assumed instantaneous in the 
analysis. H is achieved by i# ~ 0.4s which is selected as the time origin for the 
plot of Figure 3.13b. Again, sufficient time is not available to provide full agreement 
with the analytically determined similarity profile, which is shown on the figure. It is 
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Figure 3.12: h{z^ t) (mm) as a function of t (s) for six fixed locations z = Zi and 
t > 0.7s. Straight lines are fitted to the data. By interpolation, the location for 
which the slope is found to be zero is z H — 11.24mm, where H exp = 3.73 ± 0.15mm. 
The theoretical value H — 3.67mm is shown by the dashed line. 


apparent from Figure 3.13b that, after an initial overshoot, the traces are converging 
to the analytical result as time increases (see arrow on figure). In Figure 3.14, the 
time origin is shifted to t Q — 0,09s as suggested by both the inertial time scale t 0 , 
Figure 3.9, and the fit coefficient for £(£), Figure 3.8. A collapse of the data is seen 
being especially good near the tip, but a general collapse of the surface profiles for 
0.09 < t < 0.4s is also established. The similarity solution of eq 2.64 is provided for 
reference. 


3.4.2 Surface Profile, h(z,t) for data set T2 

Figures 3.15-3.18 present the data for test T2. For this test z H = 13.72mm (Figures 
3.15 and 3.16), and H exp is found to be 1.97 ± 0.08mm (Figure 3.16) compared to 
H = 1.94mm calculated by eq 2.100. The time at which the constant height condition 
is first established is found to be t H ^ 0.4s. As for test Tl, the surface traces do not 
collapse in the time available when plotted against z£ -ly/2 , Figure 3.17a. However, the 
plot of h as a function of z{t - t H )~ 1/2 , Figure 3.17b, displays nearly full agreement 
with the similarity result by the end of the test, t = 2.031s. Due to the limited 
duration of the test, ample time was not available to guarantee this observation for 
t in- The initial overshoot and convergence of the scaled profiles toward the 
theoretical solution as time increases is also observed for T2 using the time origin 
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Figure 3.13: (a) h (mm) as a function of z/t l/ 2 (mm/s 1/2 ) for nine of the ten surface 
profiles of test Tl. Trace at t = 0 is not shown, (b) As in (a) except that the 
time origin is shifted by t H -------- 0.4s. The last four traces are shown, the first trace 

is at t = 0.547s. The traces converge towards the theoretical curve as t increases as 
indicated by the arrow. 
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Figure 3.14: h (mm) as a function of z/(t — (mm/s 1 ^ 2 ) for T1 with t Q = 0.09s 

providing tip collapse. Full theory provided along with linear tip solution. The first 
trace shown is at t = 0.156s. 


in, as is indicated in Figure 3.17b. The tip location data for T2 (see Figures 3.8 and 
3.9) are best fit for the time origin t 0 = 0.03s which is employed in Figure 3.18. 


3.5 Discussion 


For cases T1 and T2, Table 3.3 lists the important time constants arising from the 
assumptions inherent in the asymptotic analysis (see Section 2.7.1). The Ohnesorge 
number Oh is also provided. The first two time scales are the viscous and slender 
body scales calculated from eqs 2.107 and 2.105 with a = 30°, 0 — 0. The third is the 
‘system’ inertial time scale which characterizes the time required for the transition 
from a gravity dominated to a surface tension dominated interface [32]. R^J 2v is the 
global viscous time scale where Ri C is the radius of the inscribed circle for the test 
cell cross-section. Comparing the magnitude of these scales, it is not surprising that 
T2 should agree more favorably with the predictions than should T1 over such short 
duration tests. Considering the low values of the time constants in the table, what 
is surprising is that T2 is not fully predicted by the theory. 3 

3 Cases T1 and T2 were selected from a family of tests conducted to explore repeatability and 
sensitivity to container surface condition. Within each family, whether the containers were initially 
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Figure 3.15: Eleven traces of the surface height, h(z,t) measured along x for T2. As 
time increases the profiles recede at the centerline x & 25mm and advance along the 
wall at x = 31.57mm, The surface profiles are shown for times t = 0, 0.039, 0.078, 
0.117, 0.156, 0.195, 0.586, 0.976, 1.367, 1.758, and 2.031s. (a) Global profile to cell 
centerline, (b) enlarged view of dashed box area on (a). The coordinates x and z are 




Figure 3.16: h(z u t) (mm) as a function of t (s) for T2. Each data set is at constant 
Zi where there is a constant increment of approximately 0.8mm. 
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Figure 3.17: (a) h (mm) as a function of z/t 1/2 (mm/s 1 / 2 ) for ten of the eleven surface 
profiles of test T2. Trace at t = 0 is not shown, (b) As in (a) except that time origin 
is shifted by tn — 0.4s. The last five traces are shown, the first trace is at t = 0.586s. 
The traces converge towards the theoretical curve as t increases, as indicated by the 
arrow. 
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Figure 3.18: h (mm) as a function of z/(t — t 0 ) 1 ^ 2 (mm/s 1 ^ 2 ) for T2 with t Q = 0.03s 
providing tip collapse. Full theory provided along with linear tip solution. The first 
trace shown is at t = 0.0391s. 


Recall that for comparison with the asymptotic analysis, z is measured relative 
to Zff and that for both T1 and T2 t H ~ 0.4s. However, fluid flows into the corner 
(z > 0) prior to tn and so there is an initial volume of fluid within z. > 0 by the 
time t = tff. This is the reason that the scaled traces of Figures 3.13b and 3.17b first 
overshoot and then converge towards the predicted behavior. As t — tn increases, the 
initial volume of fluid which was present in the corner at time tn constitutes a smaller 
portion of the total volume in the corner and the traces begin to show evidence of 
collapse about the predicted profile. The time necessary for complete collapse over 
the entire domain z > 0 is presently uncertain, particularly in light of the fact that 
t H — 0.4s while the duration of the experiment was only about 2s. 

Using t 0 as determined from the inertial reorientation time of the interface, the 
scaled traces of Figures 3.14 and 3.18 reveal collapse at the tip for t > 2 t Q and show 


dry or were purposefully pre-wetted with a thin film of the test fluid prior to the test, the mea- 
surements of the surface profiles were statistically equivalent. This agreement is not common, but 
occurs for these tests because of the natural pre-wetting of the corner prior to the drop tests despite 
the presence of gravity. This is not the case for all cells, but for tests T1 and T2 the presence of a 
thin pre-wetting film along the corner should classify the entire family as employing the ‘pre-wet’ 
surface condition. 
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Test 

Oh 

i 

H 2 / Su 

2 

10.7 /iH /a 

3 

(pRlM 1/2 

4 

Rll tv 

T1 

0.008 

0.8 

0.004 

0.12 

10.5 

T2 

0.05 

0.05 

0.05 

0.05 

0.6 


Table 3.3: The Ohnesorge number Oh and time constants, in seconds, for T1 and 
T2. 

remarkable linearity for the traces taken at times t < tH- 


3.5.1 Tip Solution Approximation 

Figures 3.14 and 3.18 show many surface profiles for which the height is nearly a 
linear function of zj(t — to) 1 / 2 . These linear profiles are reminiscent of the leading 
order tip solution identified in Section 2.6.1 where h = l — rj/y/2 and rf tip = \/2- This 
approximate solution is compared to the data for T1 and T2 on Figures 3.14 and 
3.18, respectively. The experimentally determined tip location is nearly bracketed 
by the full similarity and linear tip solutions. In fact, for test T1 which was of 1.01s 
duration, the linear tip solution agrees well with the data while the full solution does 
not. That such agreement is observed for the linear tip solution is less likely due 
to coincidence than it is to the strong influence of the solution at the tip, especially 
prior to formation of a system length scale such as H. The flow velocities are highest 
at small times. As the column extends in length, the time for information to be 
propagated downstream from z = 0 to the tip is t ~ £?. This creates a time lag 
between the origin z — z H and the tip. Consequently, though the constant height 
boundary condition is established at t = in, the tip will not respond to this condition 
until some time later. Thus, the tip region, for a time, acts independently of changes 
in conditions at z — 0, and the tip solution is better suited for predictions of the flow 
at small times. 


3,5.2 Start-up Problem 

As typified in Figure 3.18 the linear surface profiles for t < tu increase in height, 
then for t ~ t H the profiles collapse, and then give way to the predicted form of 
the similarity solution as t increases above t#. The fact that the linear profiles of 
these surfaces are reminiscent of other solutions to the governing similarity ODE (see 
Figure 2.8) makes one suspect that this ‘start-up 7 problem might also be analyzed. 

Recalling the transformation h — T a F(rj), it is surprising that for a = 1/5 a 
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Figure 3.19: h/t l/b (mm/s 1/5 ) as a function of z/t 3/5 (mm/s 3/5 ) for T2 for the first five 
traces (t <t H ) excluding the trace at t = 0. The uncollapsed trace is for t = 0.039s 
and /i/£ 1/5 = 2.62 at the origin. The full theoretical and linear tip solutions are 
provided. 


favorable collapse of the surfaces for t < t H is observed for T2 as shown in Figure 
3.19. The case of a = 1/5 is coincidentally the case of constant flow at rj — 0 solved 
in Section 2.6.3, the exact solution of which is presented on the figure. 4 A cause for 
the slight overprediction of the data by the theory could be attributed to dynamic 
contact effects, since velocities perpendicular to the contact line (and thus Ca) are 
maximum during this portion of the flow (see Section 2.7.3). The value of the flow 
rate Q is necessary to quantify such effects as will be illustrated. 

The linear tip solution for the constant flow problem also well predicts the surface 
profile during start-up as shown on Figure 3.19. Figure 3.20 plots h(zi , t) for T2 on a 
In-ln scale for the Zi nearest zh as determined by interpolation (as described in Section 
3.4.1). The slope of 1/5 is indicated as are the times for the first two surface traces 
appearing in Figure 3.18 and 3.19. The constant height time tn is also identified. 

The fact that the initial transient of the constant height solution is well governed 
by, in this case, the constant flow problem, suggests that the time constants of the 


4 The condition of a — 1/5 was not found to be true for all of the tests evaluated in this manner, 
but was found to vary somewhat with initial Bond number and fluid properties. The case of T2 is 
selected for ease of discussion. 
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Figure 3.20: h(z h t) (mm) as a function of t (s) for T2. Data taken at z* ~ z H - H , 
tn and slope a — 1/5 are indicated. 


latter situation supersede in relevance those of the former. The time constants for 
the constant flow problem are listed in Table 3.4 for tests T1 and T2. The constants 
are functions of Q and Qti ^ Qt 2 - The first constant is the viscous time scale (eq 
2.125) while the second stems from the e 2 < 1 constraint (eq 2.123). The latter is 
completely negligible in contrast to the former which may help to explain the quality 
of the agreement between the data and predictions in Figure 3. 19. For T2 it is possible 
to estimate Q using h/t 1 ^ == 2.62 (mm/s 1 / 5 ) at z = 0. This gives Q ~ 7.8 x 10 -8 
(m 3 /s) and the viscous time scale can be calculated to be ~ 1.5s. This implies that 
at z = 0, the worst case, inertia is always present during the constant flow process. 
(Recall that the constant height condition is established at approximately t ~ 0.4s 
for this test.) Inertial effects decay as ~ h 2 as the tip is approached which further 
vindicates the linear tip solution. An estimate of the capillary number using eq 2.130 
with Q for T2 shows Cad — 5.3 x 10“ 5 t“ 4 / 5 which, employing Hoffman’s correlation 
eq 2.110, produces dynamic contact angle effects on the order 6d ~ 20°, and 10° at 
times t — 0.039 and 0.4s, respectively. The impact of a change in effective contact 
angle of this magnitude is moderate and / increases 14% for 6 = = 20° from 

its value at 0 = 6 stat — 0. The increase in / (decrease in surface curvature) is not 
expected to decrease the tip velocity by more than 8% (see eq 2.114 adapted to the 
constant flow scenario). 


70 



Test 

i 

7.54(fiQ 2 / cr) 2 ' ,s i/~ 5 / 3 

2 

96.5 {n/afVQ 1 ' 2 

T1 

T2 

5 x 10 9 Q 4/3 
1 x 10 9 <> 4/3 

2.86Q 1 / 2 

30.6Q 1 ' 2 


Table 3.4: Time constants for T1 and T2, in seconds, as determined by constraints 
for the constant flow similarity solution. 

3.5.3 Further Remarks 

As discussed in Section 2.6.4 solutions to the asymptotic eq 2.73 at the tip axe linear. 
Also, in the near tip region several of the modeling assumptions ( e.g . low inertia and 
low ^-curvature) are best satisfied. In fact, the linear tip solution does a better job 
of predicting the tip location for capillary rise than does the full solution at small 
times. These findings, along with the observation that the start-up problem for the 
capillary rise experiments may be governed by the transformation h = r a F(r ] ), where 
a for T2 was found to correspond to the constant flow similarity solution (a = 1/5), 
leads one to suspect that the problem of capillary rise might be dissected into start-up 
(t < t H ) and constant height (t > t H ) regimes. An intermediate regime must provide 
matching for the two solutions for times t ~ t H . In this intermediate regime, the 
tip, acting (temporarily) independently of changes at the origin, maintains its depen- 
dence on the constant flow-like similarity solution for start-up while at the origin the 
constant height similarity profile is established which propagates downstream. An 
intermediate region is necessary to provide matching. Numerical solutions of the full 
Navier-Stokes equations would be helpful in this respect and would aid in identifying 
the regimes through which the flow transitions enroute to the t 1 ^ 2 dependence. 

Since the capillary rise problem at small times is shared between a constant flow- 
like solution during start up, where C ~ t 3 / 5 , and the constant height solution, where 
C ~ t 1 / 2 , it is no wonder that C ~ t 1 / 2 gives a reasonable representation of the tip 
location for all time! 
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Chapter 4 


Conclusion 


4.1 Summary of Results 

4.1.1 Analytical Contributions 

In this work the governing equations describing capillary flows in containers with 
interior corners are analyzed under the constraints of a slender fluid column, slight 
surface curvature along the flow direction, low inertia, and low gravity. The quasi- 
steady viscous streamwise flow is solved, averaged over the cross-sectional area, and 
using a local mass balance leads to a second order nonlinear diffusion equation for 
the height of the meniscus h(z, t) as measured from the corner axis. 

The unique velocity scale used in the nondimensionalization of the problem con- 
verts a numerical constant of 0(1) to oo used by other investigators [9] [11] [17] to an 
0(1), tightly banded function 1/8 £ F { < 1/6 for the shear free surface condition. 
This feature of the scaling clarifies the geometric influence of corner half-angle a 
and contact angle 9 and shows that, in general, contact angle effects are manifest 
primarily through the alteration of the driving pressure gradient for the flow and not 
through significant changes in the cross-flow area. The banded nature of Fi signifi- 
cantly lessens the dependence of the analysis on any numerically computed constant 
and, thus, the effects of container size, container geometry (including aspect ratio), 
and system contact angle, whether static or dynamic, may be assessed by inspection 
of the closed form results. 

Two simple problems are considered first. The problem of a steady corner flow 
is solved to yield the maximum capillary pumping rate. Next, a solution is found 
for a perturbed fluid column of infinite extent. This solution yields a time constant 
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which is useful in estimating settling times of disturbed systems. The geometric time 
constant T is shown to increase with <5 -1 as <5 — > 0 (8 — 7r/2 — a — 0) and is minimized 
for a = 30°, 0 = 0. Gravity acting perpendicular to the corner axis is added to the 
infinite column analysis and the resulting solution provides the time response and 
stability limits of the perturbed system as a function of Bond number, Boh- 

The governing equation is then transformed into a general similarity form. For 
problems with an advancing tip, h(z tip , t) = 0, an invariant transformation is em- 
ployed which, by replacing one boundary condition, permits the use of a single 
step backwards Runge-Kutta technique when numerical solution is necessary. This 
method is accurate and efficient: it not only eliminates the need for iteration, but it 
also avoids the difficulties faced by other numerical techniques (i.e. forward shooting 
and two-point boundary value methods) in resolving the solution near the tip. 

The generalized similarity equation is solved for a range of transform parameters, 
a and 6, and it is shown that a linear tip solution is valid near and at the tip for any 
choice of a. The high degree of linearity of the solutions is noted for a Z 1/5 and for 
a — 1 an analytic solution is found which produces an exact linear surface profile. 
These linear-tending solutions provide insight into the start-up flow for the capillary 
rise experiments. 

The similarity solution of Dong and Ghatzis [17] for the case of a capillary flow in 
an initially empty corner, where at time t = 0 + fluid is introduced such that the height 
of the meniscus h( 0, t) = H is fixed for all time, is re-formulated and solved. This is 
called the capillary rise problem in this study. The linear tip solution applied to this 
problem, though only appropriate in the near tip region, also produces the correct 
constant height condition at z = 0. This feature is exploited in comparisons with 
the results of drop tower experiments for this problem. The analysis is extended by 
incorporating the influence of gravity perpendicular to the corner axis. The stability 
limit of the flow to an unfavorable orientation and magnitude of gravity is found to 
agree with that determined for the infinite column solution. The constant height 
condition is calculated using the approach of de Lazzer and Langbein [29] for the 
cases of n-sided regular polygonal and rectangular cross-sectioned containers. These 
results are applied to the experiments which employ test cells of like geometry. 

New similarity solutions for a spreading drop (constant volume) and constant flow 
condition at h( 0, t) are found, the former yielding an exact analytic expression. A list 
of closed form expressions for key characteristics of the flows for the three similarity 
solutions is provided (Section 2.6.7). Observation of this list is useful for comparisons. 
For instance, the time dependence for the tip location C for each similarity solution 
is: capillary rise C ~ t 1 / 2 , spreading drop C ~ t 2 / 5 , and constant flow £ ~ t 3 / 5 . 

The constraints on each of the three similarity solutions are quantified and also 
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listed for comparison (Section 2.7). The lubrication approximation demands a slender 
fluid column, the establishment of which requires a time scale that is proportional 
to the fluid viscosity (i. On the other hand, the neglect of inertia depends on a 
viscous time scale which is inversely proportional p. The constraints of negligible 
curvature along the flow direction, gravity acting both perpendicular and parallel 
to the flow direction, dynamic contact angle effects at the moving contact line, and 
other problem-specific conditions are also considered. 


4.1.2 Experimental Contributions 

Experiments are outlined and test data is presented for the specific problem of cap- 
illary rise in transparent containers with interior corners after a step reduction in 
gravity. Test cells are employed which minimize optical distortion due to refraction. 
A 2.2s drop tower is used to access the low-gravity environment. An extensive data 
set is collected over an unexplored range of flow parameters. Measurements reveal 
repeatability and accuracy, the role of inertia and column slenderness, and the effects 
of corner half-angle, container size, container aspect ratio, and fluid properties (with 
6 = 0). The locations of the advancing tip and global meniscus at the cell centerline 
are recorded as functions of time. An automated surface digitization technique is 
employed to measure surface profiles over the entire length of the column h(z, t). It 
is shown that though some measurements correlate well with the predictions from 
the asymptotic analysis, the full surface profile measurements uncover discrepancies. 

The initial transient due to the step reduction of gravity, characteristic of most 
drop tower tests, is governed by the inertial response time of the interface, t Q . For 
the irregular containers of this study, capillary surfaces are established when t Q ~ 
0.75(pjRf c /<7) 1 / 2 , where Ri C is the radius of the inscribed circle for the container cross- 
section. This result compares with the results of Siegert et al. [32] who study the 
response of fluid interfaces in containers of rotational symmetry. The time t 0 is useful 
in providing a time origin for the capillary dominated flow. 

It is surprising that there exists a location of the interface which appears to remain 
stationary in time during the capillary rise experiments. This observation leads to 
the identification of a system origin where the constant height assumption for the 
similarity solution is applied. This solution in turn predicts the t 1 / 2 dependence for 
the flow. Specific experimental attention is paid to the establishment of this constant 
height condition and means for determining its validity and location are presented. 
The time tjj required to establish this constant height is used as the time origin for 
comparisons of the similarity solution to the experimental results. It is shown that 
for times t >,t H the scaled surface profiles first overshoot and then collapse towards 
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the similarity solution as t increases. This ‘overshoot’ behavior is attributable to 
the fact that by the time t H a significant amount of fluid has already entered the 
corner. The period t Q < t < t H is the start-up period and for a special case (test 
T2, see Section 3.5.2) it is shown that the start-up flow is also largely governed by a 
similarity solution, but for the case of constant flow. 

Prom such results it is observed that the initial transient experienced by the fluid 
is governed first by inertia t <, t Q , then by a constant flow-like regime for t Q t ^ t# , 
where h ~ t a , then by an intermediate regime which matches the tip profile (perhaps 
determined by a similarity solution as in the case of T2) to the constant height 
similarity profile extending from the origin, and finally by the constant height solution 
alone for t > t H . For both the constant height and constant flow problems, solved 
in connection with the presentation of the data, the linear tip solution was found to 
give an adequate representation of the interface (slope and tip location) using the 
time origin of t 0 . This success is due in part to the facts that the tip solution is 
the one initially established during start-up in the absence of a system length scale, 
that for a > 1/5 all similarity profiles show a high degree of linearity, that a time lag 
exists (t ~ C 2 ) between changes in /i(0,t) and the tip, and that inertia and corner- 
axis curvature are negligible in the near tip region. Since the capillary rise problem 
at small times is shared between a constant flow-like solution, where C ~ t 3/5 , and 
the constant height solution, where £ ~ t 1/2 , it is no wonder that C ~ t 1/2 gives a 
reasonable representation of the tip location for all time! The linear tip solution is 
recommended for small times. The full similarity and linear tip solutions appear to 
bracket the measured behavior of the flow. 


4.2 Recommendat ions for Future Work 

The result of the velocity scale leading to Fj ~ 0(1) significantly lessens the de- 
pendence of the analysis on any numerically computed constant. Analytical inves- 
tigations probing the added effects of dynamic wetting, free surface shear, surface 
tension gradients due to temperature or concentration gradients, and/or a partici- 
pating overlaying fluid have a firm basis from which to proceed. 

As is seen from surface traces of the developing flows during the capillary rise 
experiments, the general form of the similarity solution is matched at the tip where 
the effects of corner-axis curvature (fe 2 < 1) and inertia are negligible. Upstream 
at the apparent origin where the constant height condition is assumed applicable, 
these constraints are not readily satisfied and unpredicted corner-axis curvature re- 
mains throughout the tests. The influence of this curvature is central to the start-up 
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problem where curvatures perpendicular and parallel to the corner-axis are of equal 
order. A numerical analysis incorporating the higher order, and presently neglected, 
curvature terms in the normal stress equation would eliminate the dependence of the 
solution on the constant height condition which is fortuitously, though incidentally, 
provided by the flow and not a true boundary condition for the governing system 
of equations. A complete numerical solution to the Navier-Stokes equations would 
add insight into the impact of inertia during the initial stages of the flow as well as 
identify the transition of the flow between certain power law regimes. 

Experiments concerning spreading drops in open corners may provide an excellent 
data base for which to compare further analyses. Since for this flow scenario an exact 
analytical expression was found, the influence of first order inertia, dynamic contact 
angle, and/or corner-axis curvature may be considered in an asymptotic analysis. 
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Appendix A 
Analytical Details 


A.l Equations and Boundary Conditions 

Governing system of equations for capillary dominated flow in an interior corner. 
All quantities are dimensional See Figure 2.1 for coordinates and notation. The 
governing equations, following Davis [33], are the incompressible Navier- Stokes and 
continuity equations 


p (v t + v • Vv) = -VP + /mV • T 

(A.l) 

T = Vv + (Vv) T 

(A.2) 

V • v = 0 

(A.3) 

where v = (u,v,w). A passive overlaying fluid is assumed and body forces are 
ignored. The boundary conditions for these equations are as follows: The no-slip 
condition is applied along the walls, v = 0 on y — ±x tana. At the free surface 
S(x,y,t), the zero stress conditions are 

T • n • tx = (1 + Sf)- 1/2 ( 1 + |VS| 2 ) -1 [(1 - S 2 + S 2 Z ) (uy + v x 

- S Z (V Z + Wy)) 

+ 2Sy (tl x -Vy- S z (u z + W x ) - S z (Vy + W*))] = 

= 0 (A.4) 

T ■ n • t 2 = (1 + 5|)- 1/2 (1 + |VS| 2 )- 1/2 [(1 - Sl){u z + w x ) 

~ Sy{V z + Wy) 

+2 S z ( u x - w z ) - S y S z (u y + v x )] = 0 

(A.5) 

where ti, t 2 , and n are outward unit tangents and unit normal to S is given by 

n = (1 + IVS'I 2 ) -1 / 2 (1, —S y , —S z ) 

(A.6) 
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ti = (1 + ^)- 1/2 (l + |VS| 2 )-V 2 (5„ 1 + S 2 Z , -S y S z ) (A.7) 

t2 = (1 + Sl)~^ 2 (S z , 0, 1) (A.8) 

t 2 is restricted to the plane y = const . for convenience. The normal stress condition 
on S 

VS 

— P + Ml • n • n = 2Ha = cV • (A.9) 

where the mean surface curvature, Ti, is given by 

2H = {1 + |V5| 2 )- 3/2 (Syy( 1 + Si) + S„(l + SI) - 2 S y S z S yz ) (A. 10) 

and where 

T • n • n = 2(1 + | VS| 2 ) -1 (u x - S y ( u y + v x ) - S z (u z + w x ) 

+ S y S Z (V Z + Wy) + SyVy + S^W ^ (A.H) 

with | VS | 2 = S y + S' 2 . The two known boundary conditions on S are S = h and 
S y = 0 on y = 0, and the two assumed conditions at the contact line (y = y ma x) are 
S — ymax / tana and 

n • k = (1 + |VS| 2 )" 1//2 (sina ± S y cosa) = cos# 

The former is the interface wall intersection condition and the latter is the contact 
angle (static) condition. Here, k = (sin a, =Fcosa, 0). 


A. 2 Fi Details 


All solutions include at least terms up to second order for completeness. Recall 

(w 0 ) = —Fih z ( note Fi > 0) 

where 


Fi 


2/ /*sin<5 rSo 
Fa Jo Jy/m 


w Q dx dy 


F a is the cross sectional area function and y = y/m, m = (tan a)/hf, 8 — 7:/2—ot—0, 
Q = tt/2 — a, <j) ~ <t> — 1 is the case of 0 = 0. 


-F/j , H 2 «l, Free Surface 
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Substituting Q for a and dividing by cos 2 a , eq 2.29 becomes 


dP 0 d 2 w 0 2 n &* w ° 
+ tan { l- 


dz dx 2 


dy 2 


(A. 12) 


subject to w 0 = 0 on x — y and 


dw 0 dS 0 dw a ^2 n _ n 

___ — tan $2 — 0 

a# ay ay 

on x — S 0 - The pressure gradient term has been rescaled, P D = P 0 / cos 2 Q, and is 
0(1). Expanding w 0 and tan 2 Q in powers of Q 2 yields 


w Q = w 0o + 0 2 u; Ol + ... 


and for eq A. 12 


dP 0 d*w 0 + {n2+ , i + .^ 


h 

3‘ 


d 2 w a 


dz dx 2 

with boundary conditions w 0 = 0 on x = y and 


% 2 


= o 


(A. 13) 
(A.14) 

(A.15) 


on x = S',. By Expanding S 0 and its derivatives in powers of Q 2 , the system of eq 
2.29 and its boundary conditions may be solved to 0(Q 2 ) giving 




(A. 16) 


W ol = W 0 „ + Pc 


Oo 


Oz 


yk 

h 


(2 y 2 - x 1 - xy) + T (x 3 _ 3 x 2 y + 2 2/ 3 ) 


4- Ayk + 


yW 

h 2 


6/1 

3h\ kh\ , 


(A. 17) 


where 


h = 4> 1 - 




Note that these expansions are correct regardless of the sign of 6 as long as 8 2 1 is 

maintained. In this way positively and negatively curved surfaces may be solved for 
cases of slight curvature. No solution of this nature is possible for |<5| ~ 0(1). The 
average velocity for the case Q 2 <1, (w 0 )n, can now be determined to 0(Q 2 ) using 

. . 2 cot tt rvm [So , \ 

(w 0 )n = W Jo jy w°i) dxdy (A. 18) 


82 



Performing the integration to 0(Q 2 ) produces 


(w 0 )n = -F n y z 


(A. 19) 


where F n is given by eq 2.45. 


3. Fjj c? Q? < 1, Captive Surface 

The solution for Fn c follows similarly as that for Fn above only the boundary 
conditions applicable to eq A. 14 are w 0 = 0 on x = y and w 0 = 0 on x = S Q . The 
zeroeth and first order solutions give 


wi = Po z 


w ° ” ~it ( x2 ~~ xy ~ $° x + s° y ) 

-jy (a: 3 - xS 2 - S Q xy - xy 2 + yS 2 + y 2 S 0 ) 
(* 2 -W- Sox + S 0 y) 


(A.20) 


(A.21) 


Substitution of the above into the integral eq A. 18 yields the full expression for Fjj c 
given by eq 2.46. 


4. F hyd , F hydc , Hydraulic Diameter Approach, Free/Captive Surface 

The average velocity for the cross flow problem may also be estimated using the 
familiar hydraulic diameter approach in which 

(™o)hyd = (A. 22) 


where 

V is the wetted perimeter, 
nondimensional result 

with 


Rh.vd — 


2 A 


(A. 23) 


1 hyd ~~ jy 

Determination of A and V geometrically leads to the 


{ w o)hyd— Fhyd q 


Fhyd — ~ 


1 / F a 


8 V / sin <5 


(A. 24) 
(A.25) 
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and Fa given by eq 2.35. For the captive surface condition the wetted perimeter, V , 
is the entire perimeter of the cross flow area such that 


Fh » d 8 (/sin <$) (l + ^f) 


(A.26) 


A. 3 Similarity Details 


In this appendix the similarity parameters a and b are determined for the cases of 
(I) capillary rise, (II) spreading drop, and (III) constant flow. 

I. Capillary Rise. 


A mathematical treatment of the generalized similarity equations and applicable 
solutions is given by [25] in which the interest was nonlinear heat conduction. Specific 
details as applicable to the fluid mechanics problem are provided below. The desired 
boundary conditions for the capillary rise problem are h( 0, r) = 1 and h(£, r) = 0; 
the equation for h(z, t) is further constrained by 


-2 




(A.27) 


Though the governing ODE is only second order, three conditions are necessary due 
to the introduction of the third unknown C(r). The first condition, that the height 
h is fixed at z = 0, is a conjecture given support by the experiments described in 
Section 3.4.1. The second boundary condition is the zero height condition at the 
meniscus tip. The final condition is an integral mass balance equation which states 
that the flow rate through the cross section at z = 0 over time r is equal to the volume 
of fluid in the corner from 2 = 0 to the tip at z = C(t). By using the transforms of 
eq 2.59 and applying eq 2.61, the integral equation can be written as 


2CiCj 
3 cl H- b T 1 




(A.28) 


It should be noted that the integral mass balance condition is always satisfied by 
any choice of a and b consistent with eq 2.61. To satisfy the condition h(0, r) = 1, 
however, it is necessary to choose a — 0, b — —1/2. With the selections C\ — 1 and 
C 2 = 2" 1 / 2 eq A.28 becomes 




= 0 


(A.29) 


II. Spreading Drop 
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For a drop of constant volume the boundary conditions are h z ( 0, r) = h(C, r) = 0. 
The integral mass balance 

f h 2 dz = k (A.30) 

Jo 

states that the volume integral over 0 < z < C is a constant n for all time. Substi- 
tuting the similarity transforms of eq 2.59 into the above yields 

^l T 2a-b j nHp p2 dv = K ( A .3i) 

C 2 Jo 

By applying eq 2.61 it can be shown that this equation yields similarity when a = 
— 1/5 and b = —2/5. For the selection of C\ = (/^/5)V 5 and C 2 = (5 2 k) - 1 / 5 eq A. 31 
simplifies to 

r F 2 dr) = 1 (A. 32) 

Jo 

r] tip can be determined directly by substituting the exact solution for F into the above 
integral (or by substituting values for F + , etc.). 

III. Constant Flow 


As for the spreading drop problem, the similarity parameters a and b are deter- 
mined from the integral mass balance. For this case 



(A.33) 


which states that the constant flow rate through the cross sectional area at z = 0 at 
time r is equal to the volume of fluid in the domain 0 < z < C(t). Substituting from 
eq 2.59 yields 

C± r 2 a -b-i r up F 2 dr] 2q (A .34) 

C 2 Jo 

Applying eq 2.61 it can be shown that similarity is achieved with a = 1/5 and 
b — —3/5. For C\ ~ (5^ 2 ) 1/s and C 2 — (5 3 ^) _1//5 , the equation can be written as 



(A.35) 


A. 4 Similarity with g x 


The influence of gravity in the ^-direction can be incorporated into the analysis. In 
the dimensional problem the pressure needs only to be modified as 



a 

hTf 


+ P9 h' 
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In nondimensional terms this produces 


P Oz = {ji + B °h) 

so that 

{wo} = -Fi (l + h 2 Bo H ) h z 

Substitution of this result into eq 2.23 noting h 2 oc A and r = Fit / 2 yields eq 2.75. 
For the infinite column problem, solving for P Qz in a similar manner and by expanding 
h = h Q + e/ii and taking h 0 — const yields eq 2.58. 



Appendix B 


Correction for Optical Distortion 


An approximate correction to optical distortions attributable to mismatched indices 
of refraction between the cell material (N^) and test fluid (iYp 2 ) could be made by a 
ray trace analysis. Assuming the surface S(y, h(z , £)) is a portion of a circle in the x-y 
plane, a simple relationship may be determined between the true value of h and that 
of the apparent h measured experimentally, here labeled h exp . The sketch of Figure 
B.l is appropriate for the long working distances of these experiments. Knowing 
the slope and intercept for the line segment O-P, the conditions of its tangential 
intersection with S Q (see eq 2.32) produce the relationship 


h = h 


where 


1 -f tan 7 tan a 
erp 1 + / — / sec 7 

(B.l) 

. ( N Dl sinaA 

-sin — 77 

V No* J 

(B.2) 


Optical distortion increases with increasing a and for the capillary rise test a maxi- 
mum was experienced for the square acrylic cell, N Dl — 1.491, and the 0.65cs PDMS 
fluid, N D2 = 1.374. This worst case gives h = 0.919 h exp , approximately a 10% cor- 
rection. For the same fluid in an equilateral triangle a correction of 5% results and 
many test corrections were 3% and below. 
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